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Abstract 

One-dimensional reaction-diffusion systems are mapped through a similarity transformation 
onto integrable (and a priori non-stochastic) quantum chains. Time-dependent properties of these 
chemical models can then be found exactly. The reaction-diffusion processes related to free fermion 
systems with site-independent interactions are classified. The time-dependence of the mean particle 
density is calculated. Furthermore new integrable stochastic processes related to the Heisenberg 
XXZ chain are identified and the relaxation times for the particle density and density correlation 
for these systems are found. 
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1 Introduction 



The physics of interacting particles out of thermodynamic equilibrium is an intriguing and challenging 
problem. For important reaction-diffusion processes which display a wide variety of interesting and 
unexpected phenomena, see |]T], ^ [], |4], |5| and references therein. They are not only of experimental 
importance but also pose a theoretical problem of considerable difficulty as there is no general framework 
for their treatment, such as detailed balance for equilibrium systems. While the classic approach, 
involving phenomenological rate equations, is sufficient to the describe the time-dependent behavior of 
non-equilibrium systems in sufficiently high spatial dimensions, in low dimensions the diffusive mixing 
is ineffective and new phenomena appear. 

As an example, consider a system consisting of a single species of particles A undergoing diffusion 
and subject to the two chemical reactions A + A — > A with rate r c and A + A — > inert with rate r a . For 
this system one expects, starting from an initial particle density p , at late times an algebraic decay of 
the mean particle concentration |J 

p A {t) ~ Ar y ; t -> oo (1.1) 

In one-dimensional systems, one has V = \ independently of the branching ratio r c /r a as opposed to 
y = 1 obtained from a (mean field) rate equation. The importance of fluctuations which invalidate the 
rate equation approach is underlined by recent experimental results for the exponent y in effectively 
one- dimensional systems. In these experiments, the particles are excitons carried by long polymers and 
the effective one-dimensionality is achieved by either constraining the polymer to move in small pores 
of a convenient medium when the excitons move between the polymers [^, || or by the structure of 
the polymer itself when the excitons move along the polymer [[|. For mixed annihilation-coagulation 
(both r a ,r c ^ 0) one finds y = 0.52 — 0.59 @] and y = 0.47(3) |§, while for the pure coagulation 
case (r a = 0), y ~ 0.48 0. In the latter case, the numerical value of the exponent y is also shown 
to be independent of the initial density po Hi, which confirms exponent universality. Furthermore, in 



that particular ID system, the amplitude A is independent of po as well and it was realized |T0[ that 
this theoretical prediction is found experimentally || as well. This means that for low-dimensional 
systems with short-ranged interactions, simple kinetic equations as first investigated by Smoluchowski 



1 1] are no longer adequate to describe the strong fluctuations occurring in these systems and leading 



to so-called anomalous kinetics [|12l . More powerful techniques to include these effects must be sought. 



Theoretical results of this nature can be obtained by a variety of means, e.g. the renormalization 
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group (see e.g. improved mean- field approaches or exact calculations, see for a collection 

of reviews. It is the aim of this paper to explore some of the possibilities of obtaining exact results. 
A convenient starting point for an exact description of a stochastic process is the formulation of the 
problem in terms of a master equation. The underlying assumption of this approach, viz. the absence 
of memory effects, is at least approximately true for many real systems. A master equation (see section 
2) is a first order linear differential equation in time for the probability distribution of the stochastic 
variables describing the system. Using standard procedures, see [14" , [TB|, [T5], IT71 it may be conveniently 
written as a matrix equation and thus its solution is turned into an eigenvalue problem^. The crucial 
point is that it has been realized in recent years that for many models of interest the matrix appearing 
in this eigenvalue problem is the quantum Hamiltonian (or the transfer matrix respectively, in the case 
of discrete time dynamics) of some integrable quantum chain [[18], [T^, ^tj, ^1], |22|, |23|, |24|, |25|, |26|, ^7], ^8, 



50, 5T|, R3L 1331] . This insight has made available the tool box of integrable models (e.g. Bethe 



ansatz |35| , |3q , |37|1 ) for these interacting particle systems far from equilibrium and has led to many 
new exact results for their dynamical and stationary properties. Clearly, the models tractable by such 
means are usually rather simple, but the observed universality of many quantities of interest make these 
models both theoretically and experimentally interesting. 

Most integrable stochastic models studied so far are directly given in terms of an integrable quantum 
chain or transfer matrix. In these cases, "just by looking" at the problem one recognizes some known 
integrable system. However, it has been observed recently that one may obtain from a given stochastic 
quantum Hamiltonian by a similarity transformation the Hamiltonian of some other stochastic process, 



see 0, |3j| [10], f40| . This relates a system which is not obviously integrable to another system which is 
and in this way exact results may be obtained. At the same time it was also realized that indeed one 
may obtain by a similarity transformation a stochastic Hamiltonian from some integrable Hamiltonian 
which does not describe a stochastic process [^5], [Ltj]. This is implicit in many earlier treatments of 
reaction-diffusion processes, but has not yet been exploited systematically. It shows that many more 
stochastic systems than previously thought are actually integrable. 

In this paper, our strategy is as follows. We start from some (a priori non-stochastic) integrable 
quantum chain H. The examples we shall consider can be analyzed using Bethe ansatz and/or free 
fermion techniques. Next we perform a similarity transformation and ask under which conditions on the 
parameters of the original integrable quantum chain H and on the transformation matrix one obtains 



3 For models with discrete-time dynamics one obtains an eigenvalue problem relating the probability distribution of 
the system at time t = n + 1 to the distribution at time t = n 
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a stochastic Hamiltonian. In this way we identify several classes of new integrable stochastic models. 
Finally, we use the knowledge provided by the integrability to obtain some exact results for these 
processes, in particular, the average density at time t, and the (longest) relaxation time in the system. 

Pursuing this program in full generality appears to be a formidable task. Therefore we shall re- 
strict ourselves to systems which are (a) translationally invariant, (b) involve particles of only one kind 
with hard core repulsion preventing double occupancy of a single site and (c) which interact only via 
effective nearest neighbor interaction. We shall make at various points comments on generalizations, 
e.g. to other boundary conditions, but we shall not work out any details in these cases. Finally, we 
shall consider only similarity transformations which are translationally invariant and which preserve 



the locality of all observables. This excludes e.g. sublattice transformations or duality transfor- 
mations E2], |33| which also give rise to transformations between different stochastic Hamiltonians or 
between non-stochastic and stochastic Hamiltonians. We would like to stress that we shall consider only 
continous-time dynamics. However, all our results apply also to models with suitably chosen discrete 
time dynamics. 

The paper is organized as follows: in section 2 we review the Hamiltonian formulation for non- 
equilibrium stochastic systems and its connection with integrable quantum Hamiltonians of magnetic 
systems. We also specify the similarity transformations which will be used to link stochastic systems 
with integrable ones and comment briefly on the relationship with continuum field theory. In section 
3, we give the full classification of the stochastic systems which can be reduced to a free fermion 
Hamiltonian with site-independent interactions. In section 4, we find stochastic systems which have the 
same spectrum as a XXZ chain with real matrix elements. In section 5, we use the link with the free 
fermion case to calculate explicitly the time-dependent particle density. Section 6 gives our conclusions. 
More technical material is relegated to the appendices. In appendix A, we extended the standard 
technique of Lieb-Schultz-Mattis for the diagonalization of a hermitian free fermion Hamiltonian to the 
non-hermitian case. Appendix B gives further details for the classification of those stochastic systems 
similar to a free fermion Hamiltonian. In appendix C we recall the relationship of one of the stochastic 
systems related to free fermions to the biased voter model. Appendix D gives some further details 
for the construction of stochastic processes from the XXZ chain. Appendix E describes the details of 
the caculation of the long-time behaviour of particle densities and two-point correlators from both the 
Bethe ansatz and the truncated equations of motion. Appendix F recalls the relationship between the 
Hamiltonian spectra for periodic and free boundary conditions. 
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2 Quantum chain formulation of stochastic processes 



In this section, we review the reformulation of a non-equilibrium stochastic system defined by some 
master equation in terms of the spectral properties of an associated quantum Hamiltonian H. To be 
specific, we only consider systems defined on a chain with L sites and two allowed states per site. We 
represent the states of the system in terms of spin configurations {cr} = {o~i, a 2 , ■ ■ ■ , &l} where a = +1 
corresponds to an empty site and a = — 1 corresponds to a site occupied by a single particle. We are 
interested in the probability distribution function P(cr; t) of the configurations {a}. Our starting point 
is the master equation for the P(a; t) 



d t P(a;t) = Y: 



w(t (x)P(r; t) - w(a r)P(a; t) 



(2.1) 



where w(t — > a) are the transition rates between the configurations {r} and {a} and are assumed to 

be given from the phenomenology. In order to rewrite this as a matrix problem, we introduce a state 
vector 

|P)=£P(Mk> (2-2) 

and eq. (|2.1| ) becomes 

d t \P) = -H\P) (2.3) 
where the matrix elements of H are given by 

(a\ H |r) = -w(t -> a) if r ^ a , (<t| H \a) = ^ w(a -> r) (2.4) 

The operator P describes a stochastic process since all the elements of the columns add up to zero. 
This conservation of probability, viz. J2a P(°"i t) = 1, is equivalent to the relation 

(s\H = (2.5) 

where (s\ = J2a (^l is a left eigenvector of H with eigenvalue 0. Correspondingly, H has at least one 
right eigenvector with eigenvalue 0. Such a vector does not evolve in time and therefore corresponds to 
a steady state distribution of the system. Since in general H is not symmetric, this steady state vector 
may be highly non-trivial. Note that all this is completely general and applies to any stochastic process 
defined by a master equation. With a view on the processes that we shall study we call H a quantum 
Hamiltonian and this formulation of the master equation the quantum Hamiltonian formalism. The 
reason for this choice of language is the fact that for the processes studied below (and, in fact, many 
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other processes as well) H is the Hamiltonian of some quantum system such as the Heisenberg XX Z 
Hamiltonian. The steady state of a stochastic system corresponds in this mapping to the ground state 
of the quantum system. 

Now, it is straightforward to translate the well-known theorems [03] about the solution of the master 



equation (2.1) to the Hamiltonian formulation at hand. In particular, the real parts of the eigenvalues Ei 
of H are non-negative, tflEi > E = 0, and there are no purely imaginary eigenvalues, which excludes the 
possibility of undamped oscillations. Furthermore, if the initial distribution P ({o~}) = P(cr; 0) = (a\ Pq) 
satisfies < P ({cx}) < 1 and (s\ P ) = 1, it follows that since 

|P> = exp(-Ht) \P ) (2.6) 

the relations < (a\P) < 1 and (s\ P) = 1 hold for all times t. This guarantees the probabilistic 
interpretation. Then time-dependent averages of an observable F are given by the matrix element 

<F>(t) = y £F(a)P(a;t) = (s\F\P) = (s\Fexp(-Ht) \P ) (2.7) 

CT 

Conversely, if a matrix H satisfies the conditions 

(s\H = , (a\H\r) < for all a ^ r (2.8) 



it follows that H is the generator of a stochastic process p4| . These conditions will play a major role 
later on. We shall refer to the first as probability conservation condition and to the second as positivity 
conditions. 

In what follows, we shall be mainly interested in averages of particle numbers rij at site j and their 
correlators. These can be expressed in the quantum spin formulation in terms of the projector 

*H (* -«*)=( S !), (2J) 

and one-point and two-point functions of the rij are then expressed as 

< hj > (t) — (s\ rij \P) , < Ujfii > (t) = (s\ hj-hi \P) (2.10) 

Two basic situations are readily distinguished from the spectrum of H . If in the limit of infinite 
lattice size L — ► oo the lowest excited states have a finite gap to the ground state energy E = 0, then 
the averages (|2.10|) will approach their steady state values exponentially fast with time. On the other 



4 We stress that the structure of these matrix elements is quite distinct from expectation values in ordinary quantum 
mechanics where one always considers matrix elements between right and left states that are hermitean conjugates. 



hand, if there is in the L — > oo limit a continuous spectrum down to Eq = 0, one expects an algebraic 
approach of the correlators to the steady state. 

The recent interest in this setup comes from the integrability of quantum Hamiltonians in one 
dimension |TS|, ^ |23], |25|, |2^, As the paradigmatic example, consider the asymmetric simple 



exclusion process where particles are diffusing to the right, A + — > + A with rate 1/q and diffusing 
to the left + A — > A + with rate g. The Hamiltonian generating the time evolution of the system 
can, for free boundary conditions, be written in the form 

L-l r 



2 f , 

3=1 



(2.11) 



where a x,y,z are Pauli matrices. Through various mappings this model arises in a variety of contexts, 
it is related e.g. to the ID Kardar-Parisi-Zhang equation or to the noisy Burgers equation, see e.g. J| 
for a review. H can be analysed exactly using the Bethe ansatz []21|. It is identical to the integrable 



quantum group U q SU(2) Hamiltonian [45. 46] and gives a clear physical interpretation to the quantum 
deformation parameter q as hopping asymmetry pl[ . The quantum group symmetry may be used for 
the calculation of correlation functions ||28|| . It is well-known that the spectrum of H has a finite gap for 
q ^ 1. For q = 1 however, corresponding to symmetric diffusion, the energy levels scale with the system 
size Ei ~ L~ 2 . In contrast to conventional equilibrium systems, the spectral properties of H depend on 
the boundary conditions. For periodic boundary conditions, it can be shown that for all values of q, 
there is no gap for the infinite system and Ei ~ L~ e , where 9 varies between | and 2 depending on the 
particle density |2l], |47|, 27| . The basic mechanism behind this is explained in appendix F. 



In fact, integrable quantum Hamiltonians can be found for much more general systems. We restrict 
our attention here to models with nearest-neighbor interactions and only consider a single species of 
particles. Unfortunately, there is no standard notation for the various rates. In table |l] we give the 
possible reaction-diffusion processes together with their rates, providing at the same time a dictionary 
between the notations used by different authors. In this paper, we shall use the notation employed in 

For the time being and for purposes of illustration let us consider besides diffusion only those 
reactions which irreversibly reduce the number of particles (that is, $l,r — &l,r — v — 0). For 
temporary convenience we rescale a — > Da and define also 



n 7777" VWfn ■ V$l5r D l ;l 




3 Or equally, from the integrability of transfer matrices for discrete time dynamics |2CJ, 



diffusion to the left 
diffusion to the right 
pair annihilation 
coagulation to the right 
coagulation to the left 
death at the left 
death at the right 
decoagulation to the left 
decoagulation to the right 
birth 

creation at the right 
creation at the left 



+ A- 
A + 0- 
A + A 
A + A- 
A + A- 
A + %- 
(H + A- 
(H + A- 
A + dS- 
+ 0- 
+ 0- 
+ 0- 



> A + 
>0 + A 
^0 + 
>® + A 
*A + 
>0 + 
>0 + 
A + A 
A + A 
A + A 
+ A 
A + 



Rates defined according to reference 



D L 
D R 
2a 

1r 
1l 
S L 
S r 
Pl 
Pr 
2v 

0"L 



^32 

Oi4 
«24 
«34 

ai2 

«42 
CZ43 

&41 

a 2 i 

^31 



U>1,1 

Wl,Q 

Wl,l 
W0,1 
Wl,0 



1,0 
0,1 
0,0 
0,1 
1,0 
0,0 
0,0 

1,1 
1,1 
1,1 

0,1 
1,0 



1 10 
plO 
1 01 

1 oo 

1 01 

1 10 
plO 
1 00 
p01 
1 00 
p01 

1 11 

plO 

1 11 

p00 

1 11 
poo 

1 01 

poo 

1 10 



Table 1: Two-sites react ion- diffusion processes and their rates as defined by various authors. 



A 



\ (q + q l ) (1 + S - 7) - a , h = - {2a + 7 (q + q' 1 



(2.13) 



2 v , , j v , 2 

Note that the ratio of the left and right rates is taken to be the same for diffusion, coagulation and 
death processes. Then the quantum Hamiltonian for the system may be written in the form 



H = D(H xxz (h, A, S) + H a + H, / + H s ) 



(2.14) 



where Hxxz{h, A, 5) is the standard XX Z quantum chain, written down again for free boundary 
conditions 

L-l 



H xxz (h,A,5) = -i^^J+i + ^K+i + M^K+i- 1 



2 U 



+ hfa + - 2) - -(1 - 6) (q - g- 1 ) [a] - a 



(2.15) 



which contains the diagonal and diffusive matrix elements while the particle annihilation terms are 
contained in 

L-l 



-2a £ q 

L-l 



-2j-l 



a 3 a i+l 



h 7 = -7 q 3 ( n j a t+i + 1 1(T t n 



(2.16) 



L-l 



-5J2q 3 (q 2 (i - + <?^ + (i - %+i 



and o- ± = |(o- :B ±z(T 2/ ) are the one-particle annihilation/creation operators. 
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Remarkably, the spectrum of H is independent of the terms contained in (|2.16 ), that is spec(H) = 
spec(DHxxz{h, A, 5)) To see this, recall that the XXZ Hamiltonian conserves the number of 

particles while the reaction terms irreversibly decrease the total particle number. Thus, H can be 
written in a block diagonal form 



H 



( Ao X 5 X a 



Mi X, 



7,5 



V 



\ 



/ 



(2.17) 



where M n refers to the n-particle states and X are the reaction matrix elements. Because of the identity 



det 



A X 
B 



det A det B 



(2.18) 



it is clear that the elements of (|2.16|) do not enter into the characteristic polynomial of H . 



The phase diagram for Hxxz(h, A, 8) is known [^, 50|. For our purposes, we need the following. 
From (|2.13| ), only the portion of the phase diagram where h + A > 1 is important for us. The spectrum 
always has a finite gap when h + A > 1 , which is realized whenever 5 ^ or q ^ 1 . Then the ground 
state is a trivial ferromagnetic frozen state. The spectrum is gapless for A + h = 1, where the system 
undergoes a Pokrovsky-Talapov transition. This situation occurs for 5 = and q = 1. We have thus 
identified the cases where the model approaches the steady state exponentially (non-vanishing gap) or 
algebraically (gapless). A special point is A = where Hxxz becomes a free fermion system and thus 
in certain cases not only the spectrum may be obtained, but also correlation functions can be calculated 
explicitly (see below). 



While large classes of reaction-diffusion problems are now recognised as being integrable, see [23 



25| , |26|j for lists including systems with more than one kind of particles^], it is still far from obvious 
how to exploit the algebraic structure hidden beneath it in order to get explicit expressions for the 



desired particle number correlators for A ^ 0. A very interesting alternative was proposed in [4S 



where subsets of observables were identified for which closed equation of motion can be obtained, thus 
leading to partially integrable systems. In any case, it is important to realize that the knowledge about 
the spectrum (and, in certain cases such as the free fermion case A = 0, about correlators) does not 



6 Alternatively, this can also be seen through a similarity transformation of the form (2.19), where B is diagonal. While 
this does not change Hxxz, the rates occurring in (2.16) become arbitrary, since they contain the elements of B. Thus 
spec(iJ) must be independent of them |24[ . 

7 We remark that up to degeneracy, the spectrum of H for the processes A + ^ + A, B + $ ^ $ + B with rate 1, 
A + B, B + B — ► A with rate a and A + B — »• + with rate 2a, is the one of H X xz(a, 1 - a, 0) @. For a = 1, 
one recovers a free fermion system and closed equations of motion can be found and solved [pit , leading for example to 
PA(t) ~ p B (t) ~ iM- 
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come from the assumption that the Hamiltonian is stochastic, but from its integrability. It is therefore 
natural to ask whether one can find new integrable non-equilibrium systems from transformations of 
known integrable quantum chains. 

Specifically, we shall investigate the transformation I 



mm 



H = BHB 



-1 



B = ®B J 

j'=i 



(2.19) 



where if is a quantum Hamiltonian with known properties, H a stochastic Hamiltonian and Bj is the 
transformation matrix B acting on the site j. From now on, we shall focus on translationally invariant 
systems, i.e. we shall consider periodic boundary conditions. To study the effect of the transformation 
B, it is sufficient to consider the effect on a two-particle Hamiltonian Hjj + \, that is we write 



H 



L 



Bj,j+i 



lj+2 ■ ■ ■ 



(2.20) 



where 1^ is the unit matrix acting on the site i. For the model (|2.14j) with left-right symmetry, that is 
q = 1, we have for example 

-2a \ 



D 



( -5 -5 

1 + 5 -1 -7 

0-11 + 5 -7 

V 2(a + 7) J 



(2.21) 



Transformations of the type ( |2.19| ) have been mainly studied when both H and H represent stochastic 
systems. Besides, an equivalent way of relating stochastic systems was developed in the context of 
the Lagrangian formulation of the model (see below), which is more suitable for field theory techniques. 
In particular, these transformations allow to treat the general model ( |2.21| ), since the observables in the 



two systems are simply related |38], |39], [40], [52]. The results for the long-time behaviour of the one- 
and two-point functions for translationally invariant initial conditions are collected in table £| For the 
derivation, see appendix E. 

From the spectrum of H, we would naively expect exponential factors e~ 2k5t to be present in the 
/c-point correlator Ck, but we also see that for example algebraic prefactors are not readily predicted 
from the spectrum alone. This expectation, however, is only valid provided there are no bound states 
in the system. This is so as long the model stays sufficiently close to the case 5 = a + 7 — 1, where the 
exact spectrum of H can be found from free fermion techniques. In general, however, one can show from 

|, [Tj]] that there does exist a bound state in the two-particle sector, which 



the Bethe ansatz 53, 35 



has the energy 45 + 4 — 2 A — 2/A.0 While the results given in table |2| are valid for homogeneous initial 
8 We thank G.M. Schiitz for a useful discussion on this point. 



s 


C\ C2 




< a + 7 
> a + 7 


r l/2 t -3/2 

exp(-25t) r 1 / 2 exp(-45t) 
exp(-2cft) exp [— 4<Jt + 2(A + 1/A - 2)t] 



Table 2: Generic long time behaviour of the one-point function C\(t) = J2j < n j > (t) and the two-point 
function C^if) = J2j < n j n j+r > (t) i n the system eq. ( |2.21| ) for r finite, translation invariant initial 
conditions and finite initial particle density. The abbreviation A = 1 + <5 — a — 7 is used . 



conditions, further power law prefactors may appear for other types of initial conditions (an example 
of this will be discussed in section 5).^ 

A few comments are in order. When a, 7 are arbitrary and 5 = 0, the model can be reduced 
3£, 1C] to the cases 7 = or a = which had so far been the only ones amenable to exact treatment, 
54] , |55| . Thus the particle number correlators can be found exactly for all times. An equivalent 



sec 



transformation was subsequently given in the continuum field theory formulation of the theory . For 
5^0, the same technique allows the transformation to a system which can be interpreted as describing 
radioactively decaying particles^ which in addition undergo diffusion and either annihilation or death 
processes. Then it can be shown that in the long-time limit the equations of motion for the n-point 
functions decouple and the long-time behaviour of the correlators can be extracted DTOf. Furthermore, 
for left-right nonsymmetric rates (q 7^ 1) by the same reduction technique it can be shown that for 
$l,r — 0, the coagulation rates ^l,r can be eliminated if and only if the free-fermion condition is 
satisfied [|TIJ and the correlators are then again known exactly for all times. 

Here, we concentrate on transformations between an integrable, but not necessarily stochastic, 
quantum chain H and a stochastic Hamiltonian H . Some examples of this have already been given 
[|38| , |10|1 . We attempt to treat the transformation of type ( |2.19|) systematically starting from two types 
of an integrable quantum chainQ 

1. The most general translationally invariant Hamiltonian which be may written in terms of free 



9 The above analysis assumes that the associated amplitudes are non- vanishing. This is not always so, as can be seen in 
the calculation of the correlator < rij(t)n£ (0) > for a model in which the only non-vanishing rates are a and v. Although 
a bound state is known to be present, the associated amplitude can be shown to vanish. We thank R. B. Stinchcombe 
for pointing this out to us. 

10 The radioactive decay of particles is described in this setting by the conditions a = and 7 = S. 

11 We limit ourselves to these two types of integrable systems for pragmatic reasons. In fact, for the second type of 
integrable system we consider here we merely use the integrability of the spectrum and do not even ask whether or 
not the entire system might be integrable. A sufficient criterion to establish integrability of the full Hamilto nian H is 
provided by showing that the ifjj+i satisfy a Hecke algebra. It can be shown that for S = in the model ( 2.14 ) the 
entire non-symmetric Hamiltonian matrix and not just the spectrum-determining part is integrable |24j| . See p3| for a 
list of integrable stochastic Hamiltonians through the Hecke algebra technique. 
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fermions. We shall give in the next section the complete list of stochastic systems with spatially 
constant reaction rates which can be transformed into this form. 



2. The most general Hamiltonian with the spectrum of the XX Z chain and with real matrix elements. 
Here we shall obtain a new 3-parameter manifold (after normalization) of integrable stochastic 
processes. 

Before starting this, however, it might be useful to recall briefly how results obtained from the quan- 
tum Hamiltonian formulation of react ion- diffusion systems presented here compare with the Lagrangian 
techniques, as developed in flR O, [56]]. See |L3| for a recent review. As for the quantum Hamiltonian 



formulation, the starting point is the master equation. The creation and destruction of particles is 
conveniently described in a Fock space formalism, where the state vector becomes 

\P(t)) = $>K n 2) . . . ; t) (a\) ni (a\) n2 ■ ■ • |0) (2.22) 

cr 

where a\ creates a particle at site i and p(ni, n 2 , . . . ; t) is the probability of having n\ particles at site 
1, n 2 particles at site 2 and so on at time t. Consider for example the case when only diffusion and 
annihilation occur. Then the Hamiltonian is 

H = D J2( ajf i ~ flt i) (°* - a i) - 2a J2 ( a * 2 - flt ? a i) -> / [^(Va t )(Va) - 2a (a 2 - a f V)] d d x (2.23) 

(hi) ' " i 

where a formal continuum limit is takenP]. One then goes over to a path integral representation involving 
fields a(x,t) and a*(x,t) and an action S[a,a*]. In order to reduce the calculation of avarages to the 
usual QFT expectation values, it is customary to perform a shift a* = 1 + a. Then the action becomes 

u, o, m 

S[a,a]= J ad t a + DVaVa + Aaaa 2 + 2aa V dtd d x (2.24) 

Remarkably, since the particle number can only decrease, the diagrammatics of the associated field 
theory is greatly simplified. In particular, there are no loop corrections to the propagator and thus no 
wave function renormalization [|l^, [17], |56| and only the reaction rate a needs to be renormalized. This 
can be done by summing up the corresponding diagrams to all orders, with the result that the one-loop 
beta function is exact. From the renormalization group equations it then follows that the upper critical 
dimension is at d* = 2 and that the exponent y describing the decay of the mean particle number (|1 . 1|) 



12 In distinction to the quantum Hamiltonians written down above, one usually considers 'bosonic' operators, where 
the number of particles per site is unrestricted. While that does not appear to be an essential feature when restricting 
attention to purely destructive reactions, this may be different when considering processes which create new particles, see 

©fi- 
ll 



is given by y = min(l,d/2), in agreement with previous heuristic arguments || [j{J. In particular, it 
transpires that the whole probability distribution for fluctuations becomes universal in the late time 
regime [13| . 



What is the analogue of the similarity transformation ( |2.19j ) in the Lagrangian formulation? To see 



this, consider the process with only diffusion, annihilation and coagulation present [52|. The action 



takes the form, where V describes the interactions of the particles undergoing a reaction 

S\ = j dtdx a (a — DV 2 aj + Jdtdxdx' V(\x — x'\) [a(x)d(x') + Xa(x)] a(x)a(x') (2.25) 

such that for A = 1 only coagulation and for A = 2 only annihilation reactions occur. Following 
||52j| , consider the potentials V an for pure annihilation and coagulation, respectively, to be of the form 
V a (x) = aV(x) and Vy(x) = jV(x). Then 

2a + 7 . 
A = 2.26 

a + 7 v ' 

Thus, through the transformation a — > Xa, a — > A _1 a the coagulation-annihilation model can be reduced 

to pure coagulation. For the quantum Hamiltonian, the corresponding transformation matrix is B = 
( 1 1-A % 



A 



TO, 39 1 . In both formalisms, one arrives at the same relationship between the n-point 



correlators of the two models [|I(], [!9|, |52 

(n(xi, *i) . . . n(x n , t n )) x = (s\ a(x t , t x ) . . . a(x n , t n ) \P(t)) x 
= X-" (s\ a( Xl , t x ) . . . a(x n , t n ) \P(t)') l = X- n (n( Xl , h) . . . n(x n , t n ))[ (2.27) 

where the primes refer to the rescaling of the initial densities and correlators (since the ave rage of a 
corresponds to a particle density). 



This transformation can be extended to more general systems |60 . 

3 Stochastic systems from free fermions 

In this section we give the classification of those stochastic systems which can be obtained from a free 
fermion Hamiltonian through the transformation (|2.19|) . 

3.1 Calculation of the matrix elements 

Our starting point is the quantum chain 

L 

i=i 
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which is diagonalizable through a Jordan- Wigner transformation followed by a Bogoliubov transfor- 
mation or by a canonical transformation as described in appendix A. Here, Di 2 , f] 12 , ^1,2 are free 
parameters. In the following, we write h = hi + h 2 , g = hi — h 2 and rj 2 = r]ii] 2 . We shall consider in this 
section 77 7^ 0. The case rj = is a special case of a Hamiltonian with XXZ spectrum and will be con- 
sidered in the next section. Using (|2.20|) , we can write down the matrix describing the nearest-neighbor 
interactions 



( h-C m \ 

g-C Di 

D 2 -g-C 

V r] 2 -h-C J 



(3.2) 



and we look for a matrix H through the transformation 



L 

H = BHB- 1 , B = (g)B j (3.3) 

i=i 

and describing a stochastic system. Following the logic described in the preceding section we want 
to calculate correlation functions of the stochastic system in terms of correlators of the non-stochastic 
model which we can solve. We therefore rewrite a correlator 

(s\h Xl (ti) . ..n Xk (t k )\P ) = (s\n Xl exp (-#(ti -t 2 ))n X2 . ..n Xk exp (-Ht k ^\P ) (3.4) 

in terms of operators transformed using B, thus reducing this correlator to a matrix element computable 
through free fermion techniques. Inverting this strategy, one may first calculate matrix elements in the 
free fermion system, and then obtain the corresponding quantity for all stochastic processes related to 
this system by choosing an appropriate B. This is the strength of this approach, illustrated in section 
5 for the density p(t) as a function of time for an uncorrelated initial state with initial density po- A 
single calculation gives p(t) for all stochastic processes related to the free fermion system ( |3.1| ). 

In order to find a convenient parametrization for B we first fix the determinant of B to be 1 which 
is no loss of generality since B must be non-singular and H is independent of det B. We can then write 
B in a parametrization analogous to the Euler angles (alternatives are discussed in appendix B) 



, v a \ / cosh0 sinh0 \ / b \ , . 

\ v^ 1 ) v sinh </> c ° sh< ^ A b ~ l J 

with transformation parameters a, b, 0. Of these, we can always arrange for 6 = 1, since otherwise, 
its effect could be absorbed by redefining r\i — > b 4 r]i and r\ 2 — > b~ 4 r] 2 . So we are left with just two 
transformation parameters a and 0. Below, we shall work with Y defined by 

Y = e 20 (3.6) 
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In order to ensure that H is stochastic we shall have to satisfy the conditions ( |2.S| ). Let C m = 
J2k=i(,Hj,j+i)k,m denote the sum of the elements of the m-th column of H. Then, up to boundary 
terms in o\ L (which are absent for periodic boundary conditions) and up to an additive constant in H 



(which will be fixed later by specifying the constant C), probability conservation ( |2.8| ) is implemented 
by demanding that 

C 1 = C 4 = ^(C 2 + C 3 ) (3.7) 

These are solved by fixing 771,2 

_ (D 1 + D 2 + 2h)(l-a + Y + aY) 2 __ (D 1 + D 2 - 2h)(-l + a + Y + aYf 

Vl ~ 2(-l + a + Y + aY) 2 ' V2 ~ 2(1 - a + Y + aY) 2 { } 

The vanishing of one of the denominators is not dangerous here, because then either rji or 772 would 
have to vanish. Systems of this kind are discussed in section 4. 

For the discussion of the positivity conditions on the rates it is convenient to introduce the parity 
symmetric and parity antisymmetric combinations of reaction rates <£l,r by 

1 1 

<P = g (Pl + Vr) > <P' = 2 & L ~ Vr) ( 3 - 9 ) 

respectively, where if is taken to be a collective symbol for any member of the set 7, 8, a, D}, see 
table |I[ In order to satisfy the positivity conditions for the stochastic rates appearing in H one might 
think that all off-diagonal elements in Hjj + i have to be negative or zero. The actual conditions on these 
matrix elements are, however, weaker. To see this, it is useful to write the nearest-neighbor interactions 
in H in form of the matrix 

( a L + a R + 2V ^Sr + v -5 l - y -25 \ 

-&r + x (6R + f3 L + Dl) + x - y -D R - y 

-ol - x -D L \5 L + (3 R + D r J -x + y -y L + y 

V -2£ Pl - x -j3 R + x 25 + 7i? + 7l / 



(3.10) 

where x, y are free parameters. This matrix form corresponds to the decomposition 

= Hf}ji + x (aj +1 - aj) + y - <rf) (3.11) 

- — — 'ohvs 
where Hjj + i is the transformed matrix which in general is non-stochastic, H contains the physical 

reaction rates and is stochastic and the last two terms parametrized by x and y are divergence terms. 

For periodic boundary conditions, these extra terms cancel, while for other boundaries, they correspond 

to surface fields. 
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In this paper, we study exclusively periodic boundary conditions. Thus it is sufficient to impose the 
positivity conditons ( |2.8| ) on the parity symmetric combinations of the matrix elements of H and on 
two other combinations of transformed matrix elements in which the terms containing x, y cancel. This 
reduces the number of inequalities one has to satisfy by two. It is important to keep this distinction 
between the transformed matrix elements (which occur in Hjj + i) and the actual rates of the stochastic 
system (which occur in ifj^ s ) which have all to be positive. 

Next, we have to write down the elements of the transformed matrix i^-j+i in terms of the matrix 
elements of i^-j+i and the transformation parameters. To simplify expressions, we make use of the 
freedom to take the normalization of time scale 

D x + D 2 = 2 (3.12) 

and we shall consider the other possible case, when Di + D 2 = 0, in appendix B. Using the convention 
of ( |3.9| ), we get for the parity symmetric matrix elements 

-25 = a 2 (-l + a-Y 2 -aY 2 )-F l -N 

-7 = ± a (a 2 -l)(Y 2 -l)-F 1 -N 

-2v = a _1 (— 1 + a + Y 2 + aY 2 ) ■ F 2 ■ N 

-a = ^(a 2 -l)(Y 2 -l)-F 2 -N 
2a 

-6 = a (Y 2 -1)(1- a + Y 2 + aY 2 ) ■ F 3 ■ N 
-P = -{Y 2 -1)(-1 + a + Y 2 + aY 2 ) ■ F 3 ■ N 

-D = ^F 4 -N (3.13) 

where 

F x = h-2ah + a 2 h + AY - AaY + QhY 2 - 2a 2 hY 2 + AY 3 + AaY 3 + hY 4 + 2ahY 4 + a 2 hY 4 
F 2 = h — 2ah + a 2 h + AaY - Aa 2 Y - 2hY 2 + Qa 2 hY 2 - AaY 3 - Aa 2 Y 3 + hY 4 + 2ahY 4 + a 2 hY 4 
F 3 = h-2ah + a 2 h + 2Y - 2a 2 Y + hY 2 + 2ahY 2 + a 2 hY 2 

F 4 = h - 2ah + 2a 3 h - a 4 h + 8aY - 16a 2 Y + 8a 3 Y - hY 2 + 6ahY - 6a 3 hY 2 + a 4 hY 2 
- hY 4 - QahY 4 + 6a 3 hY 4 + a 4 hY 4 - 8aY 5 - 16a 2 F 5 - 8a 3 F 5 
+ hY 6 + 2ahY & - 2a 3 hY 6 - a 4 hY 6 
N- 1 = Y(l-a + Y + aY) 2 (-l + a + Y + aY) 2 (3.14) 
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It is remarkable that the matrix elements can be factorised in terms of only very few building blocks. 
This observation is going to be essential in the classification of the possible stochastic systems. Whether 
this factorization property has a deeper algebraic meaning is still open. Note that if we use the scaling 
\(Di + D 2 ) = —1, we recover the same expressions with Y — > —Y, which gives nothing new. The 
parity-antisymmetric matrix elements are 



i -v = 


{Dx- 


D 2 + 2g)^a(Y 2 - 


1) 


a' + x = 


Pi- 




1) 


S' + y = 


(Di- 


D 2 -2g)^a(Y 2 - 


1) 


-P'-x = 


(Di- 


D 2 2g) 1 (Y 2 
2 yJ 8aY K 


1) 


D' = 


-iPi 


-^(Y 2 + l) 





(3.15) 



In order to get a stochastic system from the matrix elements ( 3.13] , 3.15 ), we must implement the 



positivity conditions (|2.8|) on the parity symmetric rates a, 7, z/, a, 5, (3, D. These conditions completely 



determine the transformation matrix B. The physical systems thus found will have left-right symmetric 
reaction-diffusion rates. Then, in a second step and using the extra freedom coming from adding the 

^^"DiXVS ________ 

divergence terms to -ffjj+i as done in ( |3.11| ) , we shall identify free fermion systems with left-right biased 



diffusion and reaction rates. 

We now turn to the analysis of the positivity conditions (fX 



3.2 Reality of the transformation parameters a, Y 



Before embarking on checking the positivity of the symmetric rates ( |3.13| ), it is useful to break the 
problem into smaller pieces. It is a necessary condition to positivity that all matrix elements of H are 
real. This reality condition allows us to distinguish a few subcases. First, for Y 2 = 1, the transformation 



from H to H merely amounts to the diagonal transformation encountered before [5T|. So from now on 
we always assume Y 7^ 1,-1. Second, since the ratios 'j'/a' and 8' / f3' must be real, it follows from 
( FTT5D that 

a 2 real (3.16) 

and we distinguish the following. 

1. a real. Then, since 25/7 is real, it follows that 



a(-l + a-Y 2 -aY 2 ' 
(a 2 -l)(F 2 -l) 
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real (3.17) 



and it follows that Y 2 must be real. In the first case, Y itself is real, and reality of H is guaranteed 
if h, g and D\ — D 2 are real. In the second case, Y is purely imaginary and reality of the elements 
of H is assured if T) x — D 2 , h and g are all imaginary. 

2. a imaginary. We write a = iA with A real. Since 2V/a must be real, it follows from ( |3.13| ) that 

-1 + a + Y 2 + aY 2 



(a 2 -l)(F 2 -l) 

which can be rewritten as 



real (3.18) 



9 2 

Y 2 



Y 2 + 1 

^yrri = ( 3 - 19 ) 

or, letting Y 2 = P + iQ, we find that P 2 + Q 2 = 1 and thus have 

1 if a is imaginary (3.20) 

In summary, we have to consider the following cases separately 

A) a real and Y real , B) a real and Y imaginary , C) a imaginary and \Y\ = 1 (3.21) 
and we now ask when the conditions a, (3, 7, 5, V, a, D > can be simultaneously met. 
3.3 A no-go theorem 

In this section, we prove the lemma stated below. The reader not interested in the details of the proof 
can go on to the next section, where the implications will be discussed. 

Lemma: For the parity-symmetric reaction rates 5,7, v and a from eq. ( 3.19{) with the scaling h(Di + 
D 2 ) = 1, the inequalities 

a>0 , 7>0, v > , a>0 (3.22) 

cannot be simultaneously satisfied. 

Proof: We show this by examining each of the cases in ( |3.21| ) separately. For the sake of clarity 
and brevity, we are not going to go through the inequalities explicitly, but rather prefer to illustrate 
graphically what is happening. Once having understood from the graphics the result to be proved, it is 
a straightforward matter to check the relevant inequalities. The matrix elements depend on the three 
free parameters a, Y and h. We are interested in localising the regions where the rates (p (where as in 
(|3.9|), <p can stand for any of the rates 5,7, V, a) have a constant sign. These regions are bounded by 
those lines along which the (p have a simple zero or a simple pole. We now fix h and look for the curves 
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a = a(Y) such that the rates (p vanish. These curves are readily identified from the factorized form of 
the matrix elements in eq. ( 3.13| ). 
Case A: a, Y both real. 

For h fixed, the changes of sign of a are along the curves 

1 + Y 2 



a 
a± 
Y 



1 — Y 2 
h (l-Y 2 



h{l + Y 2 ) + 2Y (l ± VT^h 2 







(3.23) 



7 changes sign along the curves 



a± = 


h(l 


a = 


±1 , 


Y = 


±1 . 



h(l + Y 2 ) + 2Y (1 ± Vl - h 2 



(3.24) 



v changes sign along the curves 
l-Y 2 



a 

a± 
a 
Y 



1 + Y 2 



h-AY + 6hY 2 - AY 3 + hY 4 



-i 



l-Y 2 ) h 1+Y 



2Y(l± VT^h 2 








(3.25) 



a changes sign along the curves 

a± = [h - AY + QhY 2 - AY 3 
a = ±1 , a = 
Y = ±1 , Y = 



hY 4 



l-Y 2 ) h 1 + Y 



2Y(l± VT^h 2 



(3.26) 



Next, we have to distinguish between the cases h 2 > 1 and h 2 < 1. In the first case, the a± given above 
are complex and will not play a role in dividing the (a, Y) plane into regions where the (p are positive or 
negative, respectively. In fact, in this case the boundaries are completely independent of h. In figure 1, 
we show the excluded (shaded) regions for the four rates a, 7, V, a. Comparing them, it follows that no 
region in the (a, Y) plane remains where the four inequalities ( [3.22j ) are all satisfied. 

In the second case, the a± are all real and must be taken into account. In order to illustrate the 
mutual exclusion of the allowed regions, we show in figure 2 the excluded (shaded) regions for the four 



It 



rates for the special value h = |v3. In comparing the excluded regions for the different rates, note 
that because of the factor structure of the matrix elements, the boundary curves in the different figures 
coincide. Again, it is apparent that ( |3.22|) cannot be satisfied. 
Case B: a real and Y imaginary. 

We write Y = iy and h = ik where y, k are both real. We fix k and look for the curves in the (a, y) 
plane where the (p have simple zeroes or poles. For a, these occur at 



i + y 2 



k(l + y 2 







-1 



fc(l - y 2 ) + 2y(l ± VI + k 2 ) 



(3.27) 



7 changes sign at 



a 

y 



k(l + y 2 ) 

±1 , a = 




k(l - y 2 ) + 2y(l ± y/l + k 2 ) 



(3.28) 



v changes sign at 



i-y 2 



i + y 2 



a = 
a = 

y = o 



k - Ay - 6ky 2 + Ay 3 + ky 4 ] 1 [(1 + y 2 ) (k(l - y 2 ) - 2y(l ± y/l + k 2 ) 



(3.29) 



a changes sign at 



a = \k - Ay - 6ky 2 + Ay 3 + ky 4 ] 1 [(1 + y 2 ) (k(l - y 2 ) - 2y(l ±\/l + k 2 
a = ±1 , a = 

y = o 



(3.30) 



In analogy what was done in case A, one can show that each region in the (a, y) plane is excluded by 



at least one of the four rates and thus (|3.22|) cannot be satisfied. 
Case C: a imaginary and \Y\ = 1. 

We write a = iA and Y = cos q + i sin q with A, q real. We fix h and look for the simple zeroes or poles 
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of the rates ip. For a, these occur at 

A± = h~ l esc q ■ (l±y/T— h 2 + h cos g 
A = cot g 

For 7, they occur at 



(3.31) 



For u, they occur at 
At 



A±. = /i 1 esc g • ( 1 ± vl — h 2 + h cos g 
A = 
g = 0,±tt 



[— 3h + 4cosg — /icos2g] 1 2 sin g • f — 1 ± Vl — h 2 + h cos g 

— tang 





(3.32) 



(3.33) 



and for a, they occur at 

A± = [— 3/i + 4cosg — /icos2g] _1 2sing - 1 ± Vl — h 2 + h cos g^ 
A = 

g = 0,±tt (3.34) 

For h 2 > 1, the A±. are complex and can be ignored. For h 2 < 1, the A-t are real. To illustrate the 
mutual exclusion of the positivity conditions, we display in figure 3 the graphs for the case h 2 > 1. 
Again, it is apparent that (|3.22|) cannot be all satisfied. The case h 2 < 1 is treated similarly. 

This proves the assertion. q.e.d. 

3.4 The parity-symmetric stochastic systems 

The lemma just proven severely restricts the possibilities for the choice of the transformation parameters 
a, Y and still find a stochastic quantum Hamiltonian H. In fact, it follows that at least one of the rates 
a, 7, v, a has to vanish. But because of the explicit factorization ( |3.13| ) of the matrix elements, this 
means that two of the four rates considered will be zero. Comparing with (|3.13|) , we see that one of the 
following conditions must be implemented 



F l = , F 2 







-1 + a-Y 2 -aY* 







-1 + a + Y 2 + aY 2 







(3.35) 
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We turn to them one by one. 

1. F 1 = 0. This means that 5 = 7 = 0. We solve this for h and get 

h = — 4Y(1 — a + Y 2 + aY 2 )Ni 
iVf 1 = 1 - 2a + a 2 + QY 2 - 2a 2 Y 2 + F 4 + 2aY A + a 2 Y 4 

and the non-vanishing rates then become 

-2v = -(l + Y 2 )(l-a-Y 2 -aY 2 )^ 
a 

-5 = 2a(l-Y 2 )(l-a + Y 2 + aY 2 )N 1 
-P = 2{Y 2 - + a + Y 2 + aY 2 )N 1 

-or = -(a 2 -i)(i-y 2 )(i + y 2 )^! 

a 

-D = 2(1 + Y 2 )(-1 + a -Y 2 -aY 2 )^ 

A necessary condition for positivity are the relations 

2v D 
— > , — > 

(3 ~ 5 ~ 

Provided none of the factors in the rates vanishes, we then find 

ir 2 + i 

a Y 1 — 1 

This means that one of the following conditions must be satisfied: 



In the first case, we write a = iA with A real, i.e. Y 2 = (1 — iA)/(l + iA), and find that 



/ 4 



H 



2 + A 2 



-2A 2 



-2A 2 \ 



-2 2 + 2A 2 -2 
-2 -2 2 + 2A 2 

V o o oo/ 



does correspond to a stochastic system, which we call model I. We identify 



J\ _i_ A 2 A 2 



In the second case, we get 



£ = p = -a = 2 
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which is not stochastic. Finally, the third case Y 2 = — 1 gives, with a = %A and A real 



H 



1 + A 2 



( -2A 2 -2A 2 ^ 

2 + 2A 2 

2 + 2A 2 

V -2 -2 / 



which is again stochastic, is called model III, and we identify 



C = 1 



2A 



{A -I) 



> V2 



1 + A 2 ' 11 1 + A 2 ' ,z l + A 2 
2. F2 = 0. This means that v — a — 0. Solving for /i we find 

/i = AaY(-l + a + Y 2 + aY 2 )N 2 
N- 1 = l-2a + a 2 -2Y 2 + 6a 2 Y 2 + Y 4 + 2aY 4 + a 2 Y 4 



(3.44) 



(3.45) 



(3.46) 



The non-vanishing rates are now found as above and the discussion is completely parallel. We take 



Y — Ma - l)/(o + l) and a = iA 1 and get the stochastic quantum Hamiltonian 





H 1 

2 + A 2 



( \ 

2 + 2A 2 -2 -2 

-2 2 + 2A 2 -2 

V -2A 2 -2A 2 4 / 



(3.47) 



which is called model II and we identify 



r 1 h 9 Vl±Z A " 
G = l, /i = 2 -7— , ?7i = r?o = 7^- 

2 + A 2 ' ' ' 2 + A 2 



D 1 = D 2 



(3.48) 



This transformation of model II to free fermions had been given before ||38|| . Note that the only difference 
with model I ( |3.41| ) occurs in the sign of h, which expresses the fact that the two models I and II are 
related by a particle- hole exchange. The other possibilities (analogous to ( |3.40| )) do not give anything 
new. 

3. Y — ±y/ (a — l)/(a + 1). We write again a = iA' 1 and find in particular 

2 + A 2 



a = 4V = 2-h 



VI + A 2 



(3.49) 



This can only be made stochastic if v = a = 0. But then h = 2^1 + A 2 / (2 + A 2 ) and we are back to 
the case ( |3.47|) already found. 



4. Y = ±y(l — a)/(l + a). In the same way, we find that h = —2\/\ + A 2 / (2 + A 2 ) and we recover 
the system (|3.41|) . 
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5. a 2 = 1. We find that 8/(3 = -1. Now 8 = h(l - Y 2 )/(2Y) which only vanishes for h = or Y 2 = 1. 
In both cases, the transformation B is diagonal and we are back to the systems studied in [[HJ . 

In conclusion, we have found the complete list of parity-symmetric stochastic systems which can be 
reduced to a free fermion model with site- independent interactions through the transformation ( |2.19| ). 



These are given in ( |3.47|) , ( |3.41|) , ( |3.44|) and describe particles undergoing coagulation/decoagulation 



with the coagulation rate equal to the diffusion rate (in this case, the transformation was already known 
[BH]), the model obtained from the latter one when particles and holes are exchanged and finally the 



so-called biased voter model, see p2| . This correspondence is not immediately obvious and is explained 
in appendix C. 

3.5 Stochastic systems with a drift 

Having found the set of transformations which render the symmetric rates (p positive, we can now 
extend the identification with free fermion systems to models including a drift. 

We illustrate the procedure for the case when F\ = and then proceed to list the results. When 
F 1 = 0, we have 7 = and thus 7l = 7r = 0, since both must be positive. From (|3.15|) , this fixes 
y = (£>i- D 2 + 2g)a(l - Y 2 )/(8Y). But on the other hand, we have 

8 L + 8 R = 28 
~ ~ , ,a(Y 2 -l) 

8 L -8 R + 2y = (£>! -D 2 - 2g) K ^ } (3.50) 
Thus, solving for 8~l,r, we find 

„ a (Y 2 - 1) 

8l,r = 8 ± 1 4y ) (D 1 - D 2 ) (3.51) 

Similar considerations apply for the other rates. We can now list the stochastic systems with a drift 
which can be mapped to a free fermion system via ( |2.19|) . 

1. Model I with death/creation process ( ft.41J ) extends to the following system with the only non- 
vanishing rates 

~ 2A 2 A 2 

5 R , L = A 2 a LtR = A 2 D L>R = =F ^f^( D i ~ A,) (3.52) 

and the surface fields are given by 

1 11 A 2 

x = -(D 1 -D 2 -2g)-== , y = — (D x -D 2 + 2g) - (3.53) 

The positivity of the rates restricts D\ — D 2 to a certain interval whose boundaries depend only 
on the real parameter A. 
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2. For the coagulation/decoagulation model II ( |3.47| ) we get the following system with the non- 
vanishing rates 

~ 2A 2 A 2 

Pl,r = A 2 ^ R = A 2 Dl r = ^—^ ± —j=== (p x - D 2 ) (3.54) 

and the surface fields are given by 

x = ~ D > + 2 ^tJ=t ' y = \& - D *~ 2g) 7WTi (3 ' 55) 

For this model, the transformation to free fermions is already known p3[. 



3. Finally, the biased voter model III ( |3.44j ) is extended to the following system with the non- 
vanishing rates 

2A 2 A 

Sl,r = A 2 f3 L , R = rT -^ ± - (I* - D 2 ) (3.56) 
and the surface fields are given by 

x = -- A {D 1 -D 2 + 2g)A~ 1 , y = ~{D 1 - D 2 + 2g)A (3.57) 

In order to calculate the particle number correlators using the free fermion representation, we have 
to rewrite the particle number operator rij ( |2.9|) in the free fermion basis. This is straightforward from 
the known transformations, and thus completely determined from the parity symmetric case. First, for 
the system I ( |3.41| ), we can rewrite the parameters in the transformation matrix fl3.5| ) in the form 



A=-tan(2V>) , <j> = iij) , a = iA , 6=1 (3.58) 



and we want rij = B l rijB with n i = ^ g 1 J ' ^ e ^ nc ^ 

1 



rij = — (1 — cos 2ip a z - + sin 2tjj crj 

= ~ ( 1 ~ ~7^=< ~ , A =<%\ (3-59) 
2 V VTTa 2 ^ 3 VTTA 1 3 



Second, for the system III (|3.44 ), the transformation parameters are 



ITT 

— , a = iA , b = 1 (3.60) 



and we find 



3 2 
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Finally, for the system II ( |3.47| ), the transformation parameters are 

A = tan(2f/?) , = iij; , a = iA~ x , b = 1 



(3.62) 



and we have 

m = l 1 1 - . 1 of + . A of) (3.63) 

If one wants to calculate time correlators of the number operators in the stochastic system then one 
needs the above equalities, the state (s\B and has also to compute the action of B~ x on the initial state. 
The calculation of (s\B reduces to a single site problem since (s\ is a factorized state, (s\ = (1, 1)® L and 
also B = B® L factorizes. Later on we shall need the identity 

(1,1)5= {\fab cosh <p + v 7 a _1 6 sinh 0, -y/afe -1 sinh <p + V a _1 6 _1 cosh (/)) (3.64) 

From this one can, in particular, calculate density correlators for all four cases (see section 5). 

3.6 Summary 

To summarize our classification of the stochastic systems with site-independent interactions which 
may be transformed through a single-site transformation ( |3.3|) into a free fermion system with site- 
independent couplings, we list in table |3| the stochastic systems, as characterized by the non-vanishing 
reaction rates as defined in table [3], which can be transformed into the free fermion Hamiltonian 

H = - [ D i a J a l+i + D ^t a J+i + V (°7°i+i + a t a t+i) + h i a J + h ^ a j+i - x ] ( 3 - 65 ) 
j 

with h — hi + 1%2 and g = hi — /i2- Note that we also used the transformation parameter b in ( |3.5D and 
chosen to be 6 s = r/i/r/2 to enforce r/i = r/2 = r\ in table |3|. In table 0, we also list the values x, y for 
the surface fields which are needed to guarantee the stochastic interpretation in case of a drift, when 
D\ ^ D2, and the particle number operator fl2.9|) when transformed into the free fermion basis. 

Model I describes particles undergoing creation and diffusion with the same rate and death with an 
arbitrary rate, with and without a parity-breaking bias. Model II is obtained from this by a particle- 
vacancy exchange and describes particles undergoing coagulation and diffusion with the same rate and 
decoagulation with an arbitrary rate, with and without a drift. The transformation of this model to 



free fermions had been given before pq , p3| . Model III is the biased voter model |H| with and without 
spatial asymmetry. While the voter model is usually defined in terms of particle reactions involving 
three neighboring sites, we remind the reader in appendix C how to reformulate it in terms of two-sites 
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model 


rates / number operator 


x y 


h i] 


I 


5 R>L = A 2 d L . R = A 2 D LjR = 
A 2 [2/{A 2 + 2) t {D 1 - D 2 )/{2VA 2 + 1)] 

o n _ i 1 rr z A rr v 


-pi -p A 2 

VA 2 +1 + y/A 2 +l 


2^A 2 +1 A 2 
A 2 +2 A 2 +2 


II 


Pl,r = A 2 % tR = A 2 D LtR = 
A 2 [2/(A 2 + 2) ± (D 1 - D 2 )/(2VA 2 + 1)] 


U +^A 2 + 1 L -VA 2 +1 


2VA 2 +1 A 2 
A 2 +2 A 2 +2 


III 


$L,R = A 2 (3l, R = 

A 2 [2/(A 2 + 1) ± (D 1 - D 2 )/{2A)} 


-V+/A -V + A 


2A A 2 -l 
A 2 +l A 2 +l 


IV 


a,u,D L ,D R 
2a + 2v = D L + D R = 2, 2g = D 2 - D 1 

2n = l-a z 





a — V 2\/5^ 



Table 3: Stochastic models which can be transformed into free fermions, with the normalization (Di + 
D 2 )/2 = 1. The rates given are the only non-vanishing ones, x and y are the surface fields and 
h and r\ the parameters of the free fermion model. A is a free real parameter and the short-hand 
T>± = ±(Di — D 2 ± 2g)/4 is used. We also write the expression for the particle number operator n 
transformed into the free fermion basis n = B^ 1 nB and expressed in terms of Pauli matrices such that 
f]\ = f]2 = V in H, see eq. ( |3.65| ). 
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reactions. Finally, model IV describes annihilation and birth together with diffusion with or without 
a drift and considered in detail in ]55|, |dT| . This last model is trivial from our transformation point of 



view, because by fixing h and g the original Hamiltonian H is already a stochastic matrix itself and the 
transformation matrix B will be diagonal. 

As expected, the four models share the property that the parameters h, r\ in ( |3.65| ) satisfy h 2 +rf = 1, 
which is the well-known stochastic line in the free fermion system ( |3.65| ). These systems can therefore 
be mapped via a duality transformation |3| onto the one-dimensional kinetic Ising model with Glauber 
dynamics. However, since the transformation to free fermions of model IV is trivial and of models I, II 
and III is not, the correlators to be discussed in the free fermion language to get the particle number 
correlations are quite different. 

Another distinction becomes apparent in the presence of a drift. For models I, II, the divergence 
terms parametrized by x and y are automatically generated through the similarity transformation. 
Although for periodic boundary conditions these terms cancel, for non-periodic boundary conditions, 
they will generate surface fields. In order to find a stochastic system, it will be necessary to add into 
( |3.65|) further surface fields and see how they will transform under fl2.19|) . As a consequence, boundary 
conditions will generically affect the spectrum in an important way. This is different for model III. Here, 
for a fixed value of the parameter A, the couplings x, y of the divergence terms are proportional to each 
other and the parameter g can be chosen in ( |3.65|) such that x = y = simultaneously. Then even for 
non-periodic boundary conditions with a drift, the free fermion Hamiltonian only contains terms which 
are bilinear in the fermion operators and the techniques of appendix A for its diagonalization can be 
applied immediately. For model IV the couplings x and y always vanish. 

It is of interest to compare the list of free fermion systems of table [5] with the systems obtained 
through the approach of Peschel et al. [fpjfl . They consider a diffusive coagulation/decoagulation 
system where in addition birth and creation processes may also be present and identify certain classes 
of observables for which closed systems of equations of motion can be established and solved. The 
exitation spectrum of the corresponding quantum Hamiltonian is the one of a XY quantum chain 
in a site-dependent transverse field [i8|. Obviously, our models I, II are special cases of the systems 
studied in ||48|| , but the voter model is not included there. It remains a challenge to extend the present 
transformation technique to integrable Hamiltonians with site-dependent interactions and investigate 
which stochastic systems beyond those studied in exist (in particular, in what directions model III 
can be extended). 
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Finally, the transformation (|2.19 ) we studied is completely local. Although we derived everything 
here for one- dimensional systems with translational invariance, the transformations as they stand can 
also be applied to models which are either higher dimensional or contain non-translational invariant 
interactions, as they might arise for disordered or quasiperiodically modulated systems. Since this 
type of problem has been much studied for free fermion equilibrium models, see |)5[] for a review, the 
techniques developed for equilibrium systems can be carried over to non-equilibrium problems as well. 



4 Stochastic systems from Hamiltonian with XX Z spectrum 



In this section we give the classification of those stochastic systems which can be obtained through the 
transformation (|2.19 ) from the most general real Hamiltonian having the spectrum of an integrable 
generalized Heisenberg XX Z quantum chain. 



4.1 Determination of the transformation matrix 

Our starting point is a matrix H = ■ Hjj+i where Hj j + i is the 4x4 matrix describing the nearest- 
neighbor interactions 13 



( h u h 12 h 13 
h 22 h 23 
h 32 h 33 



h u \ 
h 2 A 



(4.1) 



V h 44 J 
and we look for a matrix H through the transformation 



H = BHB~ 1 , B = (g)Bj 

3=1 



with 



B 



bn few 
b 2 i b 22 



(4.2) 



(4.3) 



and describing a stochastic system. 

As pointed out in ]23| for a matrix of the form ( f4.1| ) the spectrum does not depend on the param- 
eters hi 2 , hi 3 , hu, h 24 , hu. Without these matrix elements, this Hamiltonian can be diagonalized using 



standard coordinate Bethe ansatz |66|, |53[. The eigenvalues of the system with L sites and periodic 
boundary conditions are 



En 



h^L 



(4.4) 



13 



One might also consider a lower triangular form for H , which physically merely amounts to a particle-hole interchange. 
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N 



E(p u ...,p N ) = E - l( 2h n ~ h 22 - h 33 ) - (h 23 + h 32 ) cosp e + i(h 23 - h 32 ) smp e ] 

e=i 

\<N<L/2 (4.5) 

N 

E(p 1 , ...,p N ) = E L - [( 2/l 44 - h 22 - h 33 ) - (h 23 + h 32 ) cosp e - i(h 23 - h 32 ) smp e ] 

e=i 

L/2 < N < L - 1 (4.6) 

E L = h M L (4.7) 



where for fixed N the pi are given by the Bethe ansatz equations 



where the ne,nj are distinct integers and the phase shifts ip PePj are given by 
'i> vm \ ( h n + ~ h 22 - h 33 ) sin ( Ei ^ £i ) 



cot 



Pe-P] 
2 



cot 



2 J (h 23 + h 32 ) cos (^) - i(h 23 - h 32 ) sin (^i) - (h n + h u - h 22 - h 33 ) cos 

(4.9) 

if 1 < iV < L/2 and 

' ip PtPj \ (^n + ^44 ~ h 22 - h 33 ) sin (gg*) 

, 2 / " [h 23 + M COS (^) + 1(^3 - /l 32 ) SU1 - (h n + h 44 - h 22 - h 33 ) COS (^) 

(4.10) 

ii L/2 < N < L — 1. From these equations we can decide how many steady states the systems has, i.e. 
how many eigenvalues 0, and whether the low-lying excitation spectrum has a non-vanishing gap or is 
gapless. 

The first step towards a classification of stochastic processes with a spectrum given by ( |4.4|) - (|4.7f ) 
consists in determining the matrix B. Since H is stochastic, Hj j + i has the left ground state 

(s\ = (1, 1)®' (4.11) 

which is a product state. This implies that also will have a left ground state as a product state and 
related to (s\ 

(s\B = (co\ = (uji,uj 2 )® 2 = (ojI,ojioj 2 ,ujioj 2 ,ojIJ (4.12) 

where uj\ t2 are free parameter. It follows that we must have 

(u\ H jd+1 = (u\ [Aj - A j+1 ) = (0, x, -x, 0) (4.13) 

where Aj is an arbitrary matrix acting on site j and x is an arbitrary constant (depending on A). Next, 
we observe from (|4.13|) that we can always arrange to have uui + u 2 — 1 . It is thus enough to keep only 



uj := lo\ as free parameter. 
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Now the implementation of eq. Q4.13 ) allows us to distinguish three cases. 



Case I: uj — 0. In this case (u\ = (0, 0, 0, 1) and by imposing ( 4.13Q we have 



whereas, from ( |4.12j ): 



B 



hdA = 



b n 1 - b 22 
-b u b 22 



Moreover, from the positivity of the spectrum we do have the following inequalities: 



h u > 



(h 23 + h 32 ) cos k + h 22 + h 33 > Vfc G [0, 2tt] 



where we have used ( |4.6|) and ( |4.10|) for N = 1. 



Case II: to — 1. The left eigenvector is (u\ = (1,0,0,0) and eq. (|4.13|) gives 



/in = 
hu + h 13 = 
h u = 



By using eq. ( }4.12|) the transformation matrix B simplifies to: 

B = 



hi -b 22 
1 - b n b 22 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



Positivity of the spectrum translates, in the same way as above, in this case into the following 
inequalities: 



flAA > 



(h 23 + h 32 ) cos k + h 22 + h 33 > V7c G [0, 2 



71 



(4.19) 



where we have used ( |4.6|) and ( |4.10|) for N — L — 1. 

Further information can be gained by considering how many steady states exist. If there is a 
single steady state, we must have 

< 6 X1 < 1. (4.20) 
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The important inequality ( |4.20| ) is seen as follows. If hu = 0, the right ground state of Hjj + i is 
|0) = (1, 0, 0, 0) = (1, 0)® 2 . This corresponds to the right ground state 0^> of Hj : j + \ through 

0) = (p,l-pf 2 = J B lJ B 2 |0) (4.21) 

where p is the mean particle density per site in the steady state and obviously < p < 1. From 
the explicit form of the matrix B, one has directly 



p = b 



ii 



(4.22) 



and (|4.20| ) follows. In addition, eq. ([4.22|) allows to calculate p once the transformation matrix is 
found. 



If there are two steady states, one of the following conditions must be satisfied: 

(h 23 + h 32 ) cos k + h 22 + h 33 = for k = 0, it 



(4.23) 



or 



h\\ = 



(4.24) 



In these cases, which we shall not treat explicitly, we did not find a similar constraint on the 
transformation parameter. 



uj 
l-u 



Case III: u / 0, 1. Since (u, 1 - u) = (1, l)K with 

11 = 

condition ( [4.13 ) can be written as 

(u\ H j:j+1 = (s\K® KHjj+iK- 1 O TV 1 



(0, x, —x, 0) 



(4.25) 



(4.26) 



In this way we have "absorbed" the parameter u in the coefficients of the Hamiltonian and 
condition fl4.13| ) simplifies to: 



l,l,l,l)H jd+1 = {0,x, -x, 0) 



(4.27) 



From eq. ( |4.26| ) we then obtain: 



h u 



-- 

hi2 + h 13 + h 22 + h 33 + h 23 + h 32 = 
hu + h 2A + h 3A + hu = 



(4.28) 
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QcllIllLlOIl 


A 


_L A WP 
1^23 "T 1I32)/ * 


o 


't 2 4 ' ' i 34 




''-22 T- 'i33 + '^23 + ''32 


p 


' 144 ^ 




g/(2P) 




Case II Case III 


5 


n/p 2 + (q + n)/p 



Table 4: Combinations of the coefficients hij that occur frequently in the text. With the exception of 
S, which is different for the cases II and III, all parameters are defined in the same way in all three 
cases. 



For B, eq. ( |4.12 ) becomes (1, 1)B = (1, 1) and we obtain 



b - ( 1 X 1 ) < 4 - 29 > 

with det-B = fen + fe 22 — 1. Positivity constraints are the same as in case II. Again, we can 
distinguish between subcases in which we have one steady steady state or two steady states 
obtaining the same additional constraints on the parameters that we obtained for case II. 

A final comment is in order. From eqs. (|4.2(j|) it is clear that H by itself also satisfies the probability 



conservation condition eq. ( |2.8|) for a stochastic system. Furthermore, the transformation ( |4.29| ) 



has exactly the same form as previously derived pi for the stochastic similarity transformation 



between two stochastic systems. Thus, taking the off-diagonal elements of H to be negative, we 
shall obtain the complete list of two-states stochastic systems which can be reduced through a 
stochastic similiarity transformation to the integrable model which just contains the rates a, 7 
and 5, corresponding to annihilation, coagulation and death processes, respectively. 

By giving the explicit forms for the transformation matrices, we have implemented the probability 
conservation condition (|2.8|) . It remains to implement the positivity conditions. As we had done before 
in the case of the free fermion integrable Hamiltonian H, we consider first the parity-symmetric case, 
which by itself will completely specify the transformation parameters fen and fe 2 2- Subsequently, we 
can then generalize to include parity-non symmetric situations, although we defer from carrying out 
this explicitly here. It will turn out that the task of analysing this set of seven coupled inequalities is 
considerably simplified by going over to new parameters, as defined in table f|. Also, it will turn out to 
be convenient in the cases II and III to replace fen and fe 22 by q% 2 as follows 
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'11 



22 



Sq 2 ; case II 

1 — fen + 5 , g 2 ; case III 



(4.30) 



4.2 Calculation of the stochastic rates 

We now give the rates of the stochastic system as obtained from the integrable chain H. 
4.2.1 Case I: integrability through equations of motion 



After performing the similarity transformation ( |4.2| ) on the matrix 

f /in h 12 hi 3 h u \ 






h 22 


h 2 3 


h 2 4 







h 3 3 




V o 








/ 



(4.31) 



with B given by ( 4.15 ) we obtain the transformed matrix H 1 . The probability conservation condition 
(|2.8|) , expressed, by eqs. (|3.7|) are here automatically satisfied and we have only to worry about the 
positivity conditions of the elements hL for all i ^ j. 

In this case, however, it turns out that the transformed elements hfj do satisfy the following pair of 
equations: 



D 1 — h 1 —h 1 —h 1 —h 1 —h 1 + h 1 + h 1 4- h 1 

LJ 1 .— f*, 34 /4 2 1 "-41 "-12 "32 ' "23 ' "43 ' "l4 
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Ml 



n3 



ii 








(4.32) 



When these conditions are both satisfied, one has a closed system of k equations for the equal-time 



fc-point correlation functions |25| . Since it can be checked immediately from the rates of the stochastic 
system whether or not the conditions eqs. ( |4.32|) are satisfied, the description either through an inte- 
grable quantum chain or through the equations of motion offer the same amount of information. While 
it might be easier to extract relaxation times from the quantum Hamiltonian, the full time-dependence 
of observables like the density is probably more easily extracted from the equations of motion. 

For practical applications it might be useful to repeat eqs. (|4.32j ) where the hij are replaced by the 
physical rates 



Pr + 1l + 2a + D R = 5 R + a R + 2u + D L 
j3 L + lR + 2a + D L = 5 L + a L + 2V + D R 

The two conditions coincide in the parity-symmetric case when they read 

5 + ct + 2? = /3 + 7 + 25 

33 



(4.33) 



(4.34) 



and we note that the value of the diffusion constant D does not enter into the integrability condition. 
4.2.2 Case II 



The matrix to be transformed is now 



H 



11 



I o h 12 -h l2 \ 

h 22 h 23 h M 

hga h 33 h u 

V h M J 



(4.35) 



and the transformation matrix B is given by (|4.18 ). Probability conservation is already implemented. 
In view of the result of the previous subsection, it is of interest to check for integrablity through the 
equations of motion. The expressions from eqs. Q4.32 ) become in this case: 



n 



D H 



- (h 12 + h 2i ) /(Sq 2 ) 
(hu - h u ) /{Sq 2 ) 



(4.36) 



and we see that the integrability from case I is only recovered if simultaneously Q = and hi 2 
(provided q 2 remains finite). 

Using the notation of table |], we can now write down the parity-symmetric rates 

P 



'31 



a 
P 
1 
5 



(gi-i)-ga (qi + R)] 

P r / 1\ 
-— 9i \Qi [Qi ~ j J ~ 9 2 ^ + R ) 

P 

—Qi m (Qi - !) - 92 {qi ~ 1 - R)] 

92 

~ (9i - l) 2 (g 2 + 1 - q x ) 
2q 2 



a 



D 



P 

-- 9i-l 

92 

P 



92 



9i - 



1 



1) - 92 (9i - 1 - R) 



9i 9i 



^) (9i - 1) - 929i (9i - 1) + ^92A 



(4.37) 



We remark that the parameters are defined such that P always emerges as a prefactor. This implies of 
course that P^O and we shall consider later what happens when P = 0. In order to get a stochastic 
system, the seven rates of eq. (|4.37|) must all be positive. We postpone the analysis of these inequalities 
to the next section and discuss first the rates for case III. 
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4.2.3 Case III 



The matrix to be transformed is 



H 



in 



( h l2 -h 12 - h 22 - h 33 - h 23 - h 32 -h 2 4 ~ hi ~ hu \ 

h 22 h 23 h 2 4 

h 32 h 33 h 34 

V flu J 



(4.38) 



and the transformation matrix B is given by (|4.29 ). Again, eqs. (3.7) are automatically satisfied with 



C\ = C 4 = and C 2 = —C 3 = h 22 + h 32 + h u . We begin by looking at the expressions D^, see ( |4.32| ), 
which for integrability through the equation of motions for the particle number correlators whould have 
to vanish and which become in this case 

D[ H = (h 23 - h 32 - h AA - hi2 - h 2A ) I (Sq 2 ) 

Dl 11 = (h 22 + 2h 32 + h 33 + h 12 - h 44 - h u ) /(Sq 2 ) (4.39) 

Now, we observe that with the parametrization of table |] and eq. (|4.30|) , the expression for rates in 
the parity-symmetric case are exactly the same as in case II and already given in eqs. Q4.37 ). Therefore 



for the analysis of the positivity, it is enough to consider one of the cases only. As a bonus, we shall 
afterwards have the freedom to work with the integrable matrices H of either case II or case III, 
depending of the applications of interest. 

We stress, however, that the coincidence of the rates for the cases II and III only holds when P ^ 0. 
When P — 0, both cases must be discussed separately. 

4.3 Analysis of the positivity conditions 

We now adress the problem of how to characterize those stochastic systems which can be transformed via 
( |2.19D to the integrable Hamiltonian H 11 or H IH . We restrict attention to the case when the elements 



of H as well as the transformation parameters qi j2 are rea/.0 We thus have to satisfy simultaneously 
the seven inequalities 

a>0 , (3>0, 7>0, 5 > , v > , a>0, D>0 (4.40) 

and the rates are given explicitly in eqs. ( f4.37|) . Since the parameter A appears only in D, we can 
always satisfy the corresponding inequality by suitably choosing A and it remains to consider the other 
six inequalities. 

14 We had already seen in the discussion of case III above that this condition might be the one needed for applications 
anyway. In addition our results from section 3, where all the elements of the free fermion Hamiltonian turned out to be 
real at the end makes it plausible that this should be the most important case. 
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While at first sight this problem looks rather difficult, it can be solved, following the ideas presented 
in section 3 after having made a few observations. From eqs. ( f4.37p , we first note that P/q 2 only appears 
as a scale factor and thus only its sign is going to be relevant for the positivity discussion. Second, 
after rescaling the rates depend only on the transformation parameters qi and q2, which characterize 
the similarity transformation B, and the parameter R, which characterizes the Hamiltonian H. We 
think of the rates as functions of these three parameters, e.g. a = a(qi, q 2 , R). We now fix R and look 
at the values of the function a in qi — q 2 space. There will be regions for which a will be positive 
and regions for which a will be negative. The boundaries of the regions G^ are characterized by the 
curves a(qi,q 2 ;R) = 0. Mapping out the boundary curves for a fixed value of R will rapidly furnish 
a geometrical idea of the shape of the region Gj°. This process is repeated for all six rates and the 
finally the desired region in qi — q 2 space where all rates are positive is the intersection 

G + = G { f n Gf n G^ ] n n G^ n G$ (4.41) 

The graphical inspection of these regions will then suggest how to satisfy in a simple way the inequalities 
(f4.40|) . Technically, the proof is much shortened by concentrating on the case with a single steady state 
which makes the inequalities ( [4.20[ ) available. Then the following distinction is sensible: i) q\ = 0, 



ii) qi = 1 and Hi) < qx < 1. (If eq. fl4.20|) does not apply, it can be shown that no further stochastic 
systems exist.) 

4.3.1 q x = 

In this case the parity-symmetric rates become 

a = 7 = 5 = 0, P = PR , + 1) 

2<? 2 

a = JL(2q 2 (l + R) + l) , D = -~PA (4.42) 
2q 2 2 

It is easy to see that for these rates, the positivity conditions can be satisfied. However, the integrability 
of this model is physically trivial, since only particle creation reactions are allowed. 

4.3.2 q x = 1 

In this case the parity-symmetric rates become 

= u = a = , a = ^(l-q 2 ) , 7 = (2q 2 (1 + R) - 1) 

5 = PR , D = --PA (4.43) 
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The discussion of positivity is completely analogous to the case above. Again, since only particle 
annihilation reactions are allowed, the integrability of the system is trivial and we just recover the 
model reviewed in section 2. 



4.3.3 < qi < 1 

We now turn to the discussion of the remaining cases. Again, the different possibilities for the signs of 
P and 92 must be treated separately. 

1. P = — 1 and g 2 > 0. From the positivity of the rates, provided that neither q\ = nor q\ = 1, we 
get the following inequalities 

a > =>- qi — q2 < 

P>0 =}► qi(qi-l)-q 2 (qi + R)>0 

7>0 =S> qi(qi-^)-q2(qi + R)>0 

S>0 =s> g 1 (g 1 -l)-g 2 (g 1 -l-i?) <0 

£ > =^ 9i - 92 - 1 < (4.44) 

5>0 (?i - - 1) -92(91 <o 

£>0 =}► 91(91 - ^)(9i-l) -92 (9i(9i-l)- A/2) >0 



From these, we observe that 



i2 < -\ (4.45) 



To see this, we note that from (3 > and 5 > we have the inequalities 91(91 — 1) > 92(91 + R) and 
91(91 — 1) < 92(91 — 1 — R), respectively. Combining these, we find that 92(91 + R) < 92(91 — 1 — -R) 
and since 92 > 0, the assertion follows. Next, we have the following 
Lemma 1: If P = —1 and 92 > and < 91 < 1, the conditions 

P > , u>0 , 5>0, D>0 (4.46) 

are sufficient for positivity. 

Proof: Since < 91 < 1, we have 91(91 — 1) < (91 — l/2){q 1 — 1). Now, a > and 5 > are 
equivalent to the inequalities (91 — l/2)(gi — 1) < 92(91 — 1 — R) and (q x — l)q 1 < 92(91 — 1 — R), 
respectively. It is now obvious that a > implies 5 > 0. On the other hand, we also have the inequality 
91(91 — 1/2) > 91(91 — 1). Now, (3 > and 7 > are equivalent to the inequalities 91(91 — 1) > 92(91 + R) 
and 91(91 — 1/2) > 92(91 + R), respectively. It follows that (3 > implies that 7 > 0. Finally, a > 
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and v > are equilavent to the inequalities gi < g 2 and q\ < g 2 + 1, respectively, which completes the 
proof. q.e.d. 
It remains to characterize G + from the remaining rates. The result is 



G+,i = { (qi,q 2 ,R) 



1 + R< qi <-R; q2 > ms *( q -i^, <^ - l '*>\ ; l) , fl - l) ; R < -\ 

\ qi + R q 1 -l-R I 2 



(4.47) 



2. P = — 1 and q 2 < 0. Now we have the following set of inequalities 

a > =^ gi — g2 > 
P>0 gi(gi-l)-g 2 (gi + i?) <0 
7>0 gi(?i-^)-g 2 (gi + i2)<0 
5>0 gi(gi-l)-g 2 (gi-l-i2)>0 

i;>0 gi-g 2 -l>0 (4.48) 

?>0 (gi-^)(gi-l)-g 2 (gi-l-i?)>0 

D>0 =^ gi ( gi _I)( gi _i)_ g2 ( gi ( gi _i)_A/2)<0 

The treatment is completely analogous to the first case. Again, from (3 and 5 we find that R < — |. 
Then we prove the 

Lemma 2: If P = —1 and q 2 < and < gi < 1, £/ie conditions 

a>0 , 7>0 , 5>0 , D>0 (4.49) 

are sufficient for positivity. 

Finally, the region G + where all rates are positive is found to be 



G+, 2 = \ (q!,q 2 ,R) 



1 + R< qi <-R-, q2 < m ^(SliSl^ t Sl(^ >qi ). R <J\ (4 . 5 o) 

V qi + R qi-l-R J 2 1 



3. P = +1 and g 2 > 0. In fact this case does not give anything new. By assumption we have 
< qx < 1. From v > it follows that 0<g 2 <gi — 1 < 0, which is impossible, because g 2 = would 
correspond to a singular transformation matrix B. 

4. P = +1 and g 2 < 0. Also here no new system arises. From a > it follows that < gi < g 2 < 
which is impossible. The only way to find stochastic systems for P positive is to consider the cases 
gi = 0, 1 which have been dealt with before. 

In conclusion, the only non-trivial parity-symmetric stochastic system which can be reduced through 
(|2.19 ) to a real integrable quantum chain with the spectrum of the XXZ Hamiltonian is the class of 
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models with rates given in eq. ( |4.37| ) and the domain of positivity 

G + = G +jl U G+, 2 (4.51) 

We remark that in the q\ — q 2 plane, the regions G +i i and G +i2 are mapped onto each other under the 
transformation (gi, q 2 ) — > (1 — qi, — q 2 ), or geometrically speaking, by point symmetry through the point 
Qi = |><?2 = 0. 

The positivity condition in ( |4.51|) is supplemented, for P = —1, but independently of the sign of q 2 , 



by the inequality 

A>2g 1 (g 1 -l) 92 ~ 9l + 1/2 (4.52) 

<?2 



Eqs. (|4.51|j4T52j ) specify completely, for the non-trivial case P = — 1, the restrictions which have to be 



put on the values of the free parameters qi, q 2 , R, A. 
4.3.4 The case P = 

We complete the analysis by stating the results in the case when P = 0. For the proofs, we refer the 
reader to appendix D. 

Case II: Here positivity requires that Q = Q = and consequently, all reaction rates vanish. The only 

processes which survive are diffusion to the right and to the left. 

Case III: If there is just a single steady state, we find the stochastic system 



( gi -hxhu -hihu ^ 

(611 - l)h u g 2 (h 23 + h 32 ) -buhu 

(6n - l)h u (h 23 + h 32 ) g 3 -buhu 

V (6n - l)h M (fen - l)hu g^ J 



(4.53) 



where the gn are determined such that the sum of each of the column vector of Hjj + % vanishes. From 
a physical point of view, however, the integrability of the model is obvious, since it describes particles 
moving diffusively and, since 7 = 5 = r_ and (3 = a = r + , decay radioactively with rate r_ and are 
created independently of their neighbors with rate r + . We remark that in this case, it is easily checked 
that the integrability conditions through the equations of motion, viz. D{ n = D 2 n = 0, are satisfied. 

4.4 Some applications 

We are now in a position to give a more physical characterization of the new integrable stochastic 
systems obtained. We shall discuss here only the class of models of which the rates are explicitly 
given in ( |4.37| ) with the domain G + ( |4 . 5 1| ) , together with ( |4.52| ) , gives the range of variables where the 
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stochastic rates are positive, and shall leave out completely all the other systems, since their behaviour 
is easily obtained without reference to the XXZ chain. 

Let us begin with a discussion of the meaning of the condition R < — |. Recalling the definition of R 
from table [| and the fact that P must be negative, we get the inequality Q > —P = Q — h 4il from which 
it follows that both /i 44 > and Q > /144 are necessarily satisfied. These are also necessary conditions 
for H to be a stochastic matrix by itself. In order to understand eq. ( |4.45| ) further, we consider the 
example of H IH in the parity-symmetric case and we also assume that H HI is itself stochastic. Then 
the matrix elements of H can be written in terms of chemical reaction rates, using the notation of 
table f| as 

h 22 = h 33 = D + 5 , h 23 = h 32 = D , h 44 = 2{a + 7) (4.54) 
We then find that Q = 25 and P = 2(a + 7 + 5) and finally 



R 



Q_ 
2P 



1 



5 



< 



(4.55) 



25-a-j ~ 2 

because Q > ^44 implies 5 > a + 7. So we see that in this particular example eq. ( ^4.45| ) simply means 
that also H is stochastic. 

Next, we ask for a characterization of integrability through similarity transformation to the XXZ 
chain solely in terms of the chemical reaction rates. In eqs. (|4.37| ), the six rates (without the diffusion 
constant D) are expressed in terms of the four parameters q±, q 2 , P and R. So we should expect to find 
two relations between the rates. Indeed, tedious calculations lead to the following relations 



-4 v 6 + 



8 5 + 4/3- 12a] v 2 + [{8a - 2(3) 5 - (p - 2d) {(3 - 6a)] v-a(p~2a 



a 



and 



+ Aa 2 V 2 - 5 2 V 2 + (28 3 + (fi - 2a) 5 2 ) u + a (-a + /?) 5 2 = 



j-4£ 3 + [47 + 45 - 12a] v 2 + 2 (7 - 2a) a 2 + [S 2 + A5a + 6 (7 - 2a) a] p] a 



(4.56) 



+ 7 2 a 2 - 5 2 a 2 + (aV 2 + 4Va + a 2 ) a 2 + {y 2 - 5 2 ) v 2 + fa 2 a + 5 3 + (7 - 2a) 5 2 ) v = 



(4.57) 



While in case I, the integrability condition could be expressed as a linear relation ([4.34j) between the 
rates, we have here two quartic relations. We note that, as in case I, the value of the diffusion constant 
D does not enter into the integrability conditions. For a given reaction-diffusion system, eqs. ( |4. 56j , fT~57| ) 
can be used to check whether or not the model is integrable through the transformation into an XXZ 
chain. 
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Since the type of integrability we consider here mainly involves the spectrum of the Hamiltonian 
H, it is particularly easy to extract the relaxation time towards the steady state and the mean particle 
density in the steady state. The relaxation time r is the inverse of the mass gap Q and thus given by 

- = Q = 2PR (4.58) 
r 

In particular, since both P and R are negative, it follows that the relaxation time r is always finite. A 
tedious calculation gives r in terms of the reaction rates 



r- 1 = -2 



(-25 - 6u + (3) a - 2v 2 + (j3 - 4?) v - 25a + (3a - 2a 2 
-2a - 4V + /?) a - 5v - 5 (5 + a) 



( 


[sv 2 + 


(-125 + 2(3 + 25] 


v + 2a 


[-2a + 


?) 


) a 


- 25 (35 + (3) 


v-25a (S + P)) 


{2av + 2v 2 + (-(3 + 4a-5) 


v - 


- a 


(-2J-6 + P) 


) 



[2aV + 2v 2 + (-(3 + 45 - 5) v - a (-25 - 5 + (3) ) 

(4.59) 

We remark that in this form there is no longer an explicit dependence on neither the coagulation rate 
7 nor the diffusion constant D. 

Finally, the mean particle density in the steady state is found from 5oo = Qi- In terms of the rates, 
it reads 



(-25 - 4V + (3) a - 5V - 5 (5 + a) 

Poo = 7 W h \ h L (4.60) 

-25 - 6V + (3 ) a - 2V 2 + ( (3 - 45) v - 5 2 - 2<55 + (3a - 25 2 



which is also not explicitly dependent on 7 or D. 
4.5 Summary 

In this section, we have found the parity-symmetric stochastic processes which are integrable through 
a similarity transformation of the type ( |2.19|) onto a real XXZ quantum chain. Up to rescaling, we find 



a non-trivial class of integrable systems, which depends on three parameters. These systems approach 
their steady steady state exponentially fast. Their rates are explicitly given in eqs. ( |Q7D and the 



constraints on the free parameters qi^, P and R are given in ( [4. 51 ). Since, working within the setting of 



case III, the integrable model is for the parity-symmetric case identical to the known system eq. ( 2.21 



we can read off from table El the time dependence of one- and two-point correlators in the t — > 00 limit. 
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5 Calculation of the density for uncorrelated initial states 
with translational invariance 



We will now illustrate the use of the similarity transformations developed in the previous sections. We 
will calculate the mean density p(t) for the various cases presented in table ^| for an uncorrelated initial 
state with density po- Since the mean density p is translation invariant, the result will not depend on 
the diffusion bias Dl — Dr. Since the bias is proportional to D\ — D 2 , it is enough to set D\ = D 2 = 1, 
which fixes the time scale as well, see eq. (|3.12|) . This calculation is greatly simplified using the free 
fermion formulation given by the Hamiltonian eq. ( [3.1D . 
The initial state is of the product form 

'*> = ®( 1 fc') J - (51) 

We wish to calculate the matrix element 

(s\ n x exp(-Ht) \p ) = (s| B n x exp(-Ht) B' 1 \p ) (5.2) 

averaged over all sites. In section 3, we had already calculated n x for the four systems of table |3] and 
related the parameters of the stochastic Hamiltonians to the parameters of the free fermion Hamiltonian 
(|3.1|) . Since the transformation ( j2.19| ) is local and \p ) is a product state, the state B^ 1 \p ) is also 
a product state. As pointed out in the previous section, the normalization of a state is arbitrary and 
therefore we can always represent \p ) = B^ 1 \p ) in the form (|5.1|) with p substituted by p . Also, 
(s\ B = (u>\, where (u>\ = ®^ =1 (oj , 1 — u>). The ratios and can be found for each separate 
case of section 3 by applying B^ 1 to |po) and B to (s\. 

The operator n x was written for models I-III in terms of the Pauli spin operators a y x and a x in 
section 3 (for model IV, trivially n x = n x ). Due to the form of (o;| we have that 

(u\ o-y = i- — - (a; I (1 - n x ) - i (u\ h x (5.3) 

(if uj 7^ 0,1). Furthermore, a x = 1 — 2h x . Therefore the calculation of the density is reduced to 
the calculation of the matrix element J2 X ( w l ™x exp(-Ht) \po). In order to do this, we write the 
Hamiltonian ( |3.1| ) in terms of spinless fermion operators by means of a Jordan- Wigner transformation, 
see appendix A. Then we express the fermion creation and annihilation operators at a given site x in 
terms of their discrete Fourier transforms to momentum space. As is well known from the properties of 
the Jordan- Wigner transformation, we have antiperiodic (periodic) boundary conditions in the sector 
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with an even (odd) number of particles, e.g. |)7, [55], |6T|. In terms of the spin language, each new 



fermion corresponds to a spin flip, the ferromagnetic state with all spins up being the fermion vacuum. 
For simplicity, we shall work in the even sector, since this simplification does not affect the results in 
the thermodynamic limit. The Hamiltonian ( |3.1| ) becomes 

H = (1 - h)L + J2 ( 2h ~ 2cos (^) ) (KK + #J-p) - 2sin (^) (V2bl p bl + »7i6p6-p) (5-4) 

where p runs over the half- integers l/2,...,L/2 — 1/2 because of the antiperiodic boundary conditions. 
The number operator is written in terms of fermion operators as J2 X ^x = J2 P bj,b p . 



Finally it can be shown | 5q , p8| that the projection of the states (lu\ and \po) over the even sector 
can be written as 

(u>r n = (0| 11(1 + J^cot (^) S P S_ P ) (5.5) 



p>0 



with f u = and 



|po>™ = AT- 1 TT (1 + /x 2 cot ( ^ ) SI^SJ) |0) (5.6) 
P>0 \ L J 

where p = p /(l — p ) and the normalization Af = ( uj \po) even . 

So we have reduced our problem to the calculation of (u\ even b p b p exp(— Ht) \po) even - Notice that 

both states (cu\ even and \po) even belong to a translation invariant subspace of the original Hilbert space 

with zero momentum. In this subspace, the operators b p b_ p , b-pbl and 1 — bU) p — St_ p S_ p are isomorphic 



to a Pauli spin algebra |6"8fl , corresponding respectively to a£, and a*. One can therefore develop 
the term exp(-Ht) using the known rules for the exponential of an expression involving Pauli spin 
matrices. First notice that 

E <^n n % h exp(-Ht) \ Po r en = 2£ />t (?£) (ur n b p L p exp(-Ht) \p )^ n (5.7) 
P P >o v L J 



since 



% + />t(^)L p ) =0 [5.1 



p 1 j uj ^ v \ L J 
which follows directly from eq. ( |5.5|) . 

Since H = J2 P >oH p , (uj\ even H = and (u\ even is a product state also in momentum space, we have 
^ | even _ g se p ara t e iy_ Therefore we can write 

E, \even ~ / ttj.\ i \even o J?2 i_ ( \ i \even it i tj±\ \ \even 

{uj\ n x exp(-Ht) \ Po ) = 2^ J u cot I — 1 {u\ b p b„ p exp(-Ht) \ Po ) 

x p >o v L J 

= 2 £ />t ( T-) ^ eV£n a v eM~H P t) \ P o) even (5.9) 
p>0 v L J 
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with Hp = 2(1 — hcosy^J) — (a p (Tp + P P <J p + 7p& p )- The constants a p , f3 p and 7 P are given by 
(Vi + V2) sm ( if) ; Pp — i (Vi ~ V2) sin ( and 7 P = (2h — 2 cos ( 77) )• ^ we define a vector a by 



Op 



the components (a p , (3 P , j p ) we have 



exp(a ■ cr p t) = cosh(| a | t) + - { <Tp sinh(| a | t) (5.10) 



a 



where cr p = (a®, cr|, cr*) and | a | = ^/a 2 + (3 p + 7 p - Substituting the definitions above and using the 
identity h 2 + 771772 = 1 one finds that | a | = 2 — 2/icos ( 77). Rexpressing a* and in terms of er+, 
cr~, one gets 

exp(a-crpt) = cosh ^(2 - 2/icos f -^r") )0 + sinh f ( 2 — 2/icos f ) ^ x 



_ cos ( 22) r/ 2 sin f ^) m sin f ^ 

1 ^ *5 + . + " . \ L 's < I .(5.11) 



Developing exp(—H p t) in ( |5.9|) with the help of (|5.11 ), using the algebra of the a matrices and the 
property of the \po) even state 

— I \ even 

0p |po) = (// tanl— tan I — J cr+ J \p ) (5.12) 
one gets after a little algebra the expression 

E, \even ~ / tt,\ i seven n v~~^ ^2 j. / TP \ / iei>en + i seven /r i o\ 

(w| n x exp(-Ht) \p ) = 2 2^ / w cot — (a; | <r+ |p ) (5.13) 

x p >o v ^ y 

4(l-hcos(^))t 

Now it can be easily shown that (uj\ even cr+ \pa) even = L ) • Substituting this expression 



in ( |5.13|) we get 

P 2 ft cot 2 < ^ 



Ei \even ~ / ttj.\ I \even r> " Jul yjyj " \ 

{u\ n x exp(-Ht) \ Po ) = 2) , , „ 2 f 2 rof 2 



f^T ( 5 - 14 ) 

p>0 ^ ' ^ Jui \ L > 



X 



, , -2 „ l ,.l-e- 4 < 1 -*'»( 2 F)>< . 



As stated above the values of f u and of p can be found from the application of the transformation 
matrix to the original states. However, one can extract f 2 from the condition (u\ even H p = 0. We obtain 
ft = > this resu h being valid for all the four cases of section 3. Separating the time-independent 



44 



part of ( |5.14| ) and using this condition, one obtains after dividing by L and taking the thermodynamic 
limit 

«o = h L dp x + pib (s) 



d gfa - (1 + e-^ 1 -^^ gin2 /g 

2n J-n P M 2j2 +tan 2^ l-/icosp Sm V2 



(5.15) 



Performing the first integral, one obtains 

1 



Pit) 



l+h 
1-h 



J_ r d Qjh - (± + h)p 2 ) e~^- hcos P^ 2 /p 

271 J-n P ^J2 +tan 2^ l-/iCOSp ^ V2 



(5.16) 



which is the central result of this section. 

When h = ±1, the generic long-time behaviour of the density is p(t) ~ p^ + At^ 1 ^ 2 , where A is 
some constant. On the other hand, when \h\ < 1, the generic long time behaviour is p(t) ~ + 
A't~ 3 / 2 exp(— t/r), where the relaxation time is 



r- 1 = 4(1 - h) 



(5.17) 



and A' is another constant. To proceed further we need the values of h, r/2, p 2 and f u for each of the 
four cases of table p|. We have 



1. For model I, we use eqns. (|3.58|J3.59| ) together with eq. ( |5.3|) to express the action of the Pauli 
spin matrices in the state (u\ in terms of n x . We obtain for the particle density 

2 _ 1 



VTTa 2 ^"' VTTa 2 

where p(t) is given by ( [5 . 1 6|) with h and 772 as a function of A were defined in section 3 and 

A 2 



(5.18) 



ft 

■J U. 



A 



P 



1 + pqVTTa^ 
1 - poVT+A? 



(5.19) 



(1 + VTTA^y 

where p is the initial particle density. Using ( |5.16|) and (|5.18|) , we get the steady state particle 
density 



1 



Pc 



(5.20) 



1 + A 2 

Model I describes particles undergoing diffusion and interacting with each other through death 
and creation processes, as defined in table [l|. Two extreme limits can immediately be understood 
intuitively. When the death processes are absent (h — > — 1 or A — > 0), the steady state density is 



45 



Poo = 1 and the approach to the steady state is algebraic. On the other hand, if only the death 
process is present (h — > or A — > oo), p^ = and the approach to the steady state is exponential. 
While in the first case, particles must diffuse to create local configurations which are favorable 
for more reactions taking place and which involves time scales of order t 1 / 2 , in the second case 
diffusion is absent and particles are just taken out of the system with a certain rate. Thus a finite 
cluster of particles will disappear within a finite time which is proportional to the cluster size and 
the death rate. 

2. For model II, where we have symmetric diffusion with coagulation and decoagulation, we repeat 
what we have done for the case above. This model had already been studied in detail before 



38], [63[ . We obtain for the density 

m = 



Pit) 



i 



VTTa 2 VTTa 2 

with p(t) given as above but with h and r} 2 given by (|3.48|) and 

r a 12 

f 2 — — 

The steady state density is 



+ 1 



(5.21) 



VTTTF + i 
A 



1 



1 + (1 - p )VT+A* 



Pc 



1 - 



1 



(5.22) 



(5.23) 



1 + A 2 

We point out that under a particle-hole transformation p — > 1 — p (and consequently p — > 1/p 
and f 2 — > 1/ we just recover the results found for model I. 

3. For model III, where we have symmetric diffusion with decoagulation and death, we repeat the 
steps of the previous cases. We obtain for the density 



2A 



A 2 



p(t) + 



1 



.4 



(5.24) 



with h and i]2 given by ( |3.45| ) and 



f 

■J u. 



A - I 
A+l 



'Ap - (1 - p ) 



(5.25) 



\Ap + (l-p )J ' 

In order to calculate the equilibrium density we have to distinguish three different cases, i.e. 
| A | < 1, | A | = 1 and | A | > 1. In the first case the steady state density is p^ = 1, this 
result being easy to understand since the rate of decoagulation (/3 = 2) is larger than the death 
rate (S = 2A 2 ). On the other hand, when | A | > 1, p^ = 0. When | A | = 1 the two processes 
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equilibrate each other and p^ = 1/2. From ( |3.45| ) we see that when A = ±1, h = ±1 and the 
decay towards the steady state is algebraic. For all other values of A the decay is exponential. 

This is in agreement with numerical studies in the biased voter model, see |62|] and references 



therein. The exponential decay of the biased voter model can be intuitively understood using the 
same argument as presented for model I. For the symmetric voter model, the algebraic decay is 



expected since this model is dual p4| to the Glauber-Ising model |TJJ at zero temperature. 



4. Finally, model IV just corresponds to model of symmetric diffusion with pair creation and anni- 
hilation. We simply have p(t) = p(t) with p{t) given by ( |5.16 ), with h = a — v and rji = 2a, 



T]2 = 2v, 2a and 2v = 2 — 2a being respectively the rates of for pair annihilation and pair creation, 
see table [3]. The constant / w = 1 and p = po/(l — p~o)- 

If we take h — ±1, then we have pure pair annihilation and pure pair creation, respectively, and 



we recover previously known results [54, 55, 38 1. The decay to the steady state of the particle 



density is algebraic and of the form t 1 ^ 2 . 

For < h < 1, we obtain an unexpected result. For p = 0, the time-dependent density had been 
calculated before j^TJ and we reproduce their result. The long-time behaviour of the density is 
generically p(t) — poo ~ t~ 1 / 2 exp(— t/r). On the other hand, if < p < 1, we get from fl5.16| ) that 
p(t) — poo ~ t~ 3//2 exp(— t/r). In addition, we have because of the particle-hole symmetry that the 
result is invariant under the transformation h — > — h, p — > So the long-time behaviour of the 
system does depend in this case on the initial condition.^ 



In summary, the main result of this section is eq. ( |5.16| ), together with eqs. (|5.18| , |5.21| , |5.24|) which 



yield the particle densities in models I, II and III. While we have limited ourselves to extract but the 
leading behaviour for large times, further terms could be readily obtained. 

6 Conclusions 

In this paper, we have addressed the problem of explicitly calculating the long-time behaviour of time- 
dependent averages of several types of reaction-diffusion systems in situations (here in low dimensions) 
where the presence of strong fluctuation effects prevents the use of simple kinetic equations. We have 
seen that, in particular in one dimension, the relationship of the master equation with the quantum 
Hamiltonians of known integrable quantum chains can be fruitfully employed. Our aim has been to 



15 



This comment also applies to models I, II and III if fj, = and 0</i<lor/i = oo and — 1 < h < 0. 
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show how to obtain, through a certain similarity transformation of the form eq. ( |2.19| ), the equivalent 
stochastic systems as characterized through their quantum Hamiltonian. Technically, the difficulty lies 
in satisfying simultaneously the set of coupled inequalities coming from the positivity of the rates which 
enter into the master equation, but from the examples presented here it should be clear how to proceed 
for more general models. 

We have obtained the following: 

1. Starting from the most general free fermion Hamiltonian with site-independent nearest-neighbor 
interactions, eq. (|3.65|) , we have found all equivalent (non-trivial) stochastic systems, as listed 
in table [3]. The temporal behaviour of the averages of all systems can be found from a single 
calculation as shown in section 5 and we have seen that a single parameter h describes the physics. 
Specifically, we have for the density pit) at large times 

f r 1 / 2 • h = l 

PV-P-~{ t-3/2exp(-t/r) 0<h<l (0) 

(provided the initial density p ^ 0) where is the steady-state density and r = — h)~ l is the 
relaxation time. Because of the particle-hole symmetry, the same form holds for the density of holes 
when — 1 < h < 0. It is known that there exist larger classes of reaction-diffusion systems which 
have the same Hamiltonian spectrum as a free fermion chain with site-dependent interactions 



54] , |4*8| . The explicit transformation generalising (|2.19| ) between these systems remains to be 
found. Another interesting extension would be the consideration of these stochastic systems on 
quasiperiodically modulated lattices. 

2. As a second example, we have started from the most general Hamiltonian which has the same 
spectrum as the XXZ chain (with site-independent interactions) and is thus treatable thorugh 
the Bethe ansatz. Up to a rescaling in the rates, we have seen in section 4 that there is a three- 
parameter family of stochastic systems equivalent to the integrable XXZ chain. The conditions for 
this type of integrability can be stated exclusively in terms of the reaction rates. Furthermore, we 
have seen that if we require parity symmetry (no left-right bias in the reaction rates) this system 
is equivalent (working with the Hamiltonian H IH in section 4) to the known reaction-diffusion 
system eq. fl2.21| ) with only irreversible reactions. The long-time behaviour of the model ( |2.21| ) 
depends on the two parameters S and A = 1 + 5 — a — 7 and is given in table |^. While for 6 = 0, 
densities and correlators decay algebraically in time, we find that for 5^0, the relaxation times 
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for multiparticle correlators depend on wether or not there exist low-lying bound states in the 
multiparticle sectors of the XXZ chain. 

This observation should remind us to be careful when extrapolating results an insights from many 
ongoing calculations, which often implicitly contain a free fermion condition (which is usually 
formulated as specifying an infinite reaction rate on the encounter of two particles), to more 
general situations, where non-trivial bound states of the quantum chain may reflect themselves in 
the values of the relaxation times. 
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Appendix A. Non-hermitian free fermion Hamiltonians 



Here we discuss how the standard techniques |57| of reducing the diagonalization of 2 L x 2 L free fermion 
Hamiltonian matrices to the diagonalization of a much smaller matrix can be extended to non-hermitian 
cases as well. For periodic boundary conditions, the complete diagonalization can be achieved by a 
sequence of Jordan-Wigner, Fourier and Bogoliubov transformations, as detailed in |55], [61], [63|]. Here 
we want a general formulation which treats all boundary conditions and space-dependent couplings on 
the same level. 

The problem of diagonalizing a non-hermitian quantum Hamiltonian such as (|3.1|) is through a 
Jordan-Wigner transformation 



m— 1 



m— 1 



(7+ = c m exp in c 



3^3 ' 



a m = exp in J2 c ) 



V c f 

u 3 u m 



(a.i; 



3=1 



reduced to diagonalizing the following quadratic form 



H 



L 

E 

n,m=l 



C n A nm C m 



(A.2) 



where A, B, D are L x L matrices. For the system ( |3.1| ) one has for periodic boundary conditions 
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h D 2 
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D = m 
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-1 

-1 

1 



V (-i) 





1 







(A.4) 

where M = Y%=\ c ln c m is the fermion number operator and the c m satisfy the usual anticommutation 
relations. For more general nearest-neighbor interactions and/or boundary conditions, the precise form 
of the matrices A, B, D will be affected but not the diagonalization technique to be described. 
Following ||67|| , we seek new fermionic variables £ q which render H diagonal 

. v l) (A ' 5) 
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and which are defined through the following canonical transformation 

& = £ ( a V C L + a _1 M gm C m ) , £ q = ( aU *qm c] m + U^qm ™) (A.6) 

m m 

where u, v are L x L matrices and a is a free parameter. For a = 1 we recover the canonical transfor- 
mation for hermitian H p7| . We stress that £j~ 7^ £t, the hermitian conjugate of £ 9 , when a/l. But 
all what we need are the anticommutation relations 

= {C.ff} = {e 9 ,^} = (A.7) 

which imply that 

^ ] (^qm^pm ^qm^prn^) ^9iP ; ^ ] i^qm^pm ^ qm^pm) (A. 8) 

We can thus use and £ 9 as particle creation and annihilation operators to generate the entire spectrum 
of H from the one-particle energies A q . Now, from (|A.5[) we must have [£ 9 , ] = A g ^ g and this leads 
to 

A.qQUq m ^ ] ^QUq n A mn OL Vq n B nm ^ , AgCl V qm ^ ] (^OlUq n D nm -\- OL Vq n A nm ^ (A. 9) 

n n 

From the structure of the Hamiltonian (|3.1|) , it is apparent that B = — (772/771) -D. We now choose 
a = (^2/ 7 7i) 1//4 - Introducing a vector notation (w g ) n = u qn , we then have 

A g (u q *,v q *) = (uq\v q *) ^ ^ J (A.10) 

which reduces the problem to the diagonalization of a 2L x 2L matrix. 

A further simplification can be achieved when A is symmetric, A = A T (which in the context of the 
model ( |3.1|) means D\ = D 2 ). We then find, analogously to |67| 



K (K + v*) = (u q * +v*)M , M = (A + a 2 D) {A - a 2 D) (A. 11) 

and have reduced the problem to the diagonalization of an L x L matrix. In particular, since M. is 
hermitian, it follows in this case that all A q are real and can be taken to be positive. The hermitian 
case treated in [H7H is recovered for a = 1. 



Appendix B. Some remarks on the transformation ( |ST3| ) 

We add some technical remarks on the transformation ( |3.3| ) used to find stochastic systems from the 
free fermion Hamiltonian Q3.1]) . 
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B.l The case D x + D 2 = 



When D\ + D 2 = 0, the analysis of section 3 has to be repeated. First, probability conservation (|2.8| ) 
is implemented by fixing i]i )2 

^±h Ma ,Y)r, ^.n^( _\- + a a + + y + ; a r ) 2 (B1 > 

The parity symmetric matrix elements are 

-2a = -a 2 h(l - a + Y 2 + aY 2 ) -G x -M 

-2D = a^hi-l + a + Y 2 + aY 2 ) ■ G 2 ■ M 

-7 = hia(a 2 - 1)(Y 2 - 1) • d ■ M 

-6 = ha(Y 2 - 1)(1 - a + Y 2 + aY 2 ) ■ G 3 - M 

-P = -h(Y 2 + a + Y 2 + aY 2 )-G 3 - M 

-a = —h(a 2 -l)(Y 2 -l)-G 2 -M 

—D = --h(a 2 - 1)(Y 2 - l) 2 ■ G 2 ■ M (B.2) 
2 

where 

Gi = 1 - 2a + a 2 + QY 2 - 2a 2 Y 2 + Y 4 + 2aY 4 + a 2 Y A 

G 2 = 1 - 2a + a 2 - 2Y 2 + 6a 2 Y 2 + F 4 + 2aY A + a 2 F 4 

G 3 = 1 - 2a + a 2 + Y 2 + 2aY 2 + a 2 Y 2 

M" 1 = Y(l-a + Y + aY) 2 (-l + a + Y + aY) 2 (B.3) 

and we see that h only enters into the matrix elements as a constant factor and thus plays no role 
whatsoever. The parity nonsymmetric matrix elements are 

-2 _ 1 
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= a- 1 (D 1 + g) 


~i 
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= a{D x + g) — 




= a-\D 1 -g) 
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2 _ 1 
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(B.4) 

2Y V ; 
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Second, we apply the reality criterion. Since both 5' / (3' and YA 7 ' must be real, it follows that a 2 is 
real. In fact, we find that the reality conditions (|3.21|) also hold here. Third, one localizes the curves 



in (a, Y) space where the parity symmetric rates change sign. Mapping out the corresponding areas, 
we find as before that the positivity conditions fl2.8|) for the rates cannot simultaneously be satisfied. 
Finally, one might try to use the factorized structure of the matrix elements in ( |B.2| ) to put selectively 
some rates to zero. Carrying out the calculation, we find that no stochastic systems occur, since there 
will be always pairs of parity-symmetric rates which are not positive multiples of each other. 
In conclusion, the case D\ + D 2 = does not lead to any stochastic system. 

B.2 Alternative forms of the transformation matrix 

Here we briefly discuss a few alternative transformation matrices. Begin with the transformation matrix 

B = ( I /-> ) < B - 5 > 

where x and y are free parameters. It is enough here to consider the unbiased case Di = D 2 = 1. 
Probability conservation is then implemented through g = and 

Vl = (h + 1)(1 + xy) 2 /x 4 , ?7 2 = (l-h)x 4 /(l + xy) 2 (B.6) 

For the rates, we find in particular 

a = -p= (h - l)xy / (I + xy) 2 (B.7) 

and, from the positivity conditions, we must have either xy = or h = 1. In the first case, the 
transformation B is singular for x = and trivial for y = 0. In the second case, f]2 = and we find 



H + l 



{ -2(1 + 2xy) \ 

1 -1 2xy 

-1 1 2xy 

V 2 J 



(B.8 

which is indeed stochastic for — | < xy < 0. We recover the annihilation-coagulation model 110, |39 



discussed in section 2. As discussed in section 3, we can include surface fields to allow for biased 
reaction-diffusion rates. Alternatively, one might try 

B = ( I T ) < R9 > 

The only stochastic form is ( [B.8| ) with y — > 1/y. Considering the inverses of the above transformation 
matrices corresponds to the particle-vacancy exchange A <-> 0. 
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Finally, we consider ( |3.1|) with D x = D 2 = D and t/i = r] 2 = T). The transformation matrix is 

b = i : 1 : ij i (b.io) 



611 612 

&21 ^22 



Let C m = Y$e=i{Hj,j+i)k,m- A necessary condition for probability conservation is C 2 — C3. This implies 



(611 + 621) (621 + 622)23 







(B.n; 



detE 

The case g = was discussed in section 3. If bu + 621 — 0, the two other conditions from probability 
conservation require that 



fc_ Hd = * _ Hd =/l+ nr 



:b.i2) 



This implies that either 6f 2 = 622 or that 77 = 0. In the first case, bi 2 = —b 22 renders B singular and 
from 6 12 = b 22) it follows 77 = anyway. So we take the second case. From 5 and 7, it then follows that 
D = as well. The resulting Hamiltonian is 



H 



( -2hV 2h 2 U 2h l U \ 

2h 2 V -2h 1 V-2h 2 U 2h 1 U 

2h 1 V -2h 1 U-2h 2 V 2h 2 U 

2h 1 V 2h 2 V -2hU J 



(B.13) 



where 



U 



>V2 



V 



>22 



(B.14) 



bl2 + ^22 b\ 2 + &22 

This is stochastic provided hi, h 2 < and 612 • b 22 > 0. Physically, the system is trivial, since it consists 
of a set of non-interacting particles which are spontaneously created and decay spontaneously. Similarly, 
if from ( |B.11| ) we have 621 + b 22 = 0, we get 77 = D = and 



H — — 



( -2hV 2h 2 U 2h 1 U ^ 

2h 2 V -2h 1 V-2h 2 U 2h 1 U 

2h 1 V -2h 1 U-2h 2 V 2h 2 U 

2h 1 V 2h 2 V -2hU J 



;b.is) 



with 



U 



bu 



V 



611 + 621 611 + 621 

which is stochastic if hi, h 2 > and bu and 621 have the same sign. The physics is again trivial. 



;b.i6) 



Appendix C. Quantum chain formulation for the biased voter 
model 



Conventionally, the biased voter model |62[ is defined as a Markov process where the sites on the lattice 
are either in state or state 1. If site j is in state 0, it changes to state 1 at a rate Ari where r\ is the 
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number of nearest neighbors of the site j which are in state 1. Similarly, sites in state 1 change to state 
at rate r Q . Thus the elementary processes involve three sites rather than two as it is in the quantum 
Hamiltonians we consider here. Now consider the stochastic quantum Hamiltonian H = J2j with 
two-body matrix 

/ -1 -1 \ 



3+1 



-1 -1 

A+l 

A+ 1 

V -A -A J 

From this we construct the three-sites contribution 



(CI) 



H j-U,j+1 = 1 ® H 3,j+1 + H 3~ 



1,3 



(C.2) 



which completely describes the transitions of the central site j. The transitions of the central site occur 
with the rates 



0A0 -> 000 2 

®AA -> UA 1 

AA$ -> AU 1 

AAA — > AtyA 



A$A^AAA 2A 

00,4 -> 0yL4 A 

A00 -> AA0 A 

000 -> 0A0 



(C.3) 



Now, if one identifies A with the state 1 and with the state 0, one reproduces for the central site j 
the rules of the biased voter model as specified above. Carrying out the sum over all lattice sites yields 
the same rules for all sites j. 

Appendix D. The case P = for the XXZ chain 

We give here the details of the calculations leading to the results stated in section 4. The cases II and 
III must be considered separately. 



D.l Case II 

The following notation is used. We write b\ = bu and 62 — ^22 and use the parameters Q and Q from 
table (|. Then the parity-symmetric rates are 
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a 



ft (61 - If 
2 b 2 



(d.i; 



D = — - 



2b 2 
1 

26 2 



2ft(6 1 --)(6 1 -l) + 6 2 Q 
2ft(6 1 -~)(6 1 -l) + 6 2 (/i 2 3 + M 



Now, we can always arrange to have ft positive. To see this, recall that the transformation matrix B 
viewed as function of 6 2 satisfies 

B(b 2 ) = B(-b 2 ) (J (D.2) 

Therefore a change of sign in b 2 changes the sign of the matrix elements h\ 2 , h 2 ^, hu- Since h\ 2 does not 
occur in eqs. ( p.l|) , we can always absorb the sign of ft into b 2 . Next, from the positivity of a and v 
one has the inequalities Qb\b 2 > and ft(6i — 1)62 > 0. The following two cases are to be distinguished: 

1. bx > 1 and b 2 > 0. From > we get 61(61 - 1) > Q/(2ft)6 2 or 61 > 1 + (Q/(2fi)(6a/&i) =: 1 +W. 
On the other hand, from 7 > we have 61(61 - |) < Q/(2ft)6 2 or 1 < 61 < \ + (Q/2ft)(6 2 /&i) = | +U. 
This implies that W > | > 0. Now, taking the two inequalities together, we have 1+W<6i<|+W 
which is impossible. The only escape is to set ft = Q = 0. 

2. 61 < and 6 2 < 0. From the positivity of 7 and a we get 2ft(6i — a;)(6i — 1/2) =p b 2 Q < 0, and 
x = 0, 1, respectively. Adding these two relations gives (61 — l) 2 < 1/2, which would imply > 61 = 1/2, 
which is impossible. Again the only escape is to set ft = Q = 0. 

But if ft = Q — 0, all rates with the exception of D vanish, which means that the only possibility 
is pure diffusion. 

D.2 Case III 

Using the same notation as above for case II, we find the parity-symmetric rates 



- 2a 
-2u 
-27 
-25 
-2/9 
-2d 
-2D 



6 2 (-6i(ft + Q))/(6i + 6 2 -l) 
(l-6i) 2 (-6i(ft + Q)+ft + Q) 

61 (2(ft + Q)b\ - 61 (ft + 2Q) - b 2 Q + Q) /(&! + 62-I) 

-61 (2(ft + Q)b\ - 6i(2ft + Q) + 6 2 g - Q) /{h + 6 2 - 1) 

(1 - 61) (2(ft + Q)b\ - 61 (2ft + 3Q) - 6 2 g + Q) /(61 + b 2 -l) 



(D.3) 



= (61 - 1) (2(ft + Q)bj - 61 (3ft + 2Q) + 6 2 Q + ft) /(6 X + 6 2 - 1 
= 26?(ft + Q) - 6 2 (3ft + 3Q) + 61 (ft + Q + h 23 + h 32 ) 
+ (6 2 -l)(/> 23 + M]/(6i + 6 2 -l) 
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For a single steady state we have < b% < 1 and we get the inequalities (Q + Q)b 1 < from 5 > and 
(Q + — 1) < from v > 0. This implies that /i 44 = Q = —ft- Inserting into the rates, one sees 

that the remaining positiviy conditions become 

(&! + b 2 - 1) > , (h 23 + h 32 ) < (D.4) 

and one recovers the stochastic Hamiltonian given in section 4. 

Appendix E. The long time behaviour of correlation functions 

In this appendix we give a brief discussion of the results presented in Table 2, namely we show how the 
long time behaviour of the one and two-point correlation functions may be obtained by a truncation of 
the equations of motion for these quantities which becomes exact at long times [ID| . It is easy to see 



that the long time behaviour of the one and two-point equal time correlation functions is determined 
by the lowest lying state in the one and two-particle sectors respectively, which is non-orthogonal to 
the initial state, provided that there is not a state in a sector with a larger number of particles that has 
a lower energy. We have seen in section 2 that the spectrum of Hamiltonian ( [2.21| ) is identical to the 
spectrum of Hxxz(h, A, 5). So the study of the spectrum of the one and two-particle sector of Hxxz 
allows one to obtain directly the time decay of these two correlation functions. In particular it can be 
shown, and we will do this below, that there is a line of bounds states in the two particle sector for any 
A^O. The presence of those states has to be taken into account in the discussion of the behaviour of 
the two point correlation functions. The energy of lowest state in the one-particle sector is, as given by 
( f4.5| ), equal to E = 26 (pi =0). So the decay of the one-particle correlation function will be determined 
by this quantity. The proof of the existence of a line of bound states in the two-particle sector is also 
straightforward. By definition a bound state between two particles, 1 and 2 say, exists when their 
momenta acquire an imaginary part. We see from equation ( |4.8| ) that such an imaginary part must be 
included in the phase shift ip PlP2 and since we are interested in the thermodynamic limit we therefore 
put ip PlP2 = i(f>L. Hence p\ = + i(f) and p 2 = ^rn 2 — icfi. Substituting these two equations in ( |4.9| ) 
with ip PlP2 of the above form one obtains 

(1 + S - a - 7 ) sin[ I {ni - n 2 ) + i0] ^ tiLft ^ ^ 



cos[ |(ni + n 2 ) ] — (1 + S — a — 7) cos[f (ni — n 2 ) + \ 2 

for L large. The minus (plus) sign applies when <j) is positive (negative), and we have substituted hn, 



h 22 , ^33, h 4 4, h 2 3 and h 32 by their values as given in (|2.21| ). We shall from now on, and as in section 
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2, use the abbreviation A = 1 + 5 — a — 7. We can, after some trivial manipulations write ( |E.1| ) as 
cos[j(rii + n 2 ) ] = Ae~\^e ±l i ( - ni ~ ri2 \ From this equation it follows that ri\ — n 2 = kn where k is an 
integer. It also follows that if | A | < 1 then the maximum value cos[ j{n 1 + n 2 ) ] can achieve is | A | < 1. 
So there is no bound state with zero total momentum within this range of A. The energy of the bound 



state can be calculated from ( }4.5|) after the appropriate substitutions and the use of the cosine addition 
formula and is equal to 

9 . r tt 1 

(E.2) 



E b = 4(1 + 8) - 2A - — cos 2 



-(m + n 2 ) 



Now one can ask when this energy will be lower than the lowest lying state in the continuum namely 
the state with energy E = 45 (pi = p 2 = 0). This will happen if 4 — 2(A + ^ cos 2 j {rii + n 2 ) ) < 
subjected to the restriction cos[ ^(ni+n 2 ) } = (— l) fc A e - '^'. From these equations one obtains A > 1 and 
E = 45 — 2(A — 2 + A^ 1 ) for the energy of the lowest lying state in the line of bound states (n\ = —n 2 ). 
So for this range of the parameter the long time behaviour of two point correlation functions will be 
determined by the energy of this state. For the other values of A the decay is determined by the state 
in the continuum with lowest energy (E = 45). These are the results presented in Table 2. 

Let us now show how these results arise in the context of the equations of motion for the correlation 



functions, following |10| . First let us consider the one-point correlation function. The equation of motion 
for the expectation value is obtained by using the Heisenberg equation of motion ^ (s\hj(t) = 



s| [H, nj(t)] where H is given by (|2.21|) and we are taking the average in the initial state in the 



way indicated by Q2.7|) . Using the properties of (s\ one obtains 

Y<nj(t)> = < n j+1 (t) > + < nj-x(t) > — 2(1 + 5)< hj{t) > 

+ (5 - 7 - 2a) ( < nj(t)n j+1 (t) > + < h^^hjit) > ) (E.3) 

where we are using the abbreviation < n,j(t) > = (s \ n,j(t) \^). Now one knows that the energy of the 
two particle lowest lying state (either the state with pi = p 2 = for A < 1 or the bound state for 
A > 1) is higher than the energy of the lowest lying state of the one particle sector (which has an 
energy equal to E = 25 (p\ = 0)). So the last term of the r.h.s of ( |E.3| ) has a faster decay than the 
remaining ones since the decay of this term is determined by a higher value of the energy. Hence we 
can drop it at large times fllQf . Assuming a translationally invariant initial state with < rij(0) > = po, 
one finds for the average particle density at t large < fij{t) > = poe~ 2St , which is the result presented in 
Table 2. The equation for the two point correlation function is obtained in a similar way. Taking the 
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expectation value of the Heisenberg equation of motion for fij(t)hi(t) (I > j) one obtains 

^-< hj{t)ni{t) > = < n j+1 {t)ni{t) > + < n^ x {t)ni{t) > + < %(t)n ; _i(t) > 



dt 

+ < hj{t)h m (t) > - 4(1 + 5)< hj(t)hi(t) > 

+ ( 5 - 7 - 2a) ( < nj(t)n j+ i(t)ni(t) > + < nj^fjhjfyn^t) > 

+ < n j (t)n l (t)n l+1 (t) > + < nj(t)ni-i(t)hi(t) > ) (E.4) 

for I ^ j + 1 and 

-^< hj(t)h j+1 {t) > = < fij-x^n^t) > + < nj(t)n j+2 (t) > - 2(1 + 5 + 7 + a)< nj(t)n j+1 (t) > 

+ ( 5 - 7 - 2a) ( < hj(t)h j+1 (t)h j+2 (t) > + < n i _i(i)fi i (t)n i+1 (t) >) (E.5) 

for / = j + 1. Now one proceeds as before. It is easy to see that there is no bound state involving three 
particles in the three particle sector. This is so because then all the momenta of the three particles would 
have an imaginary part and the energy of the state would be complex. This is impossible given that 
Hxxz-i which determines the spectrum, is hermitian. So the energy of the lowest lying state in the three 
particle sector is higher than the energy of the lowest lying state in the two particle sector. Therefore the 
three particle correlators decay faster then the two-particle correlators and can at large times be dropped 
from equations ( E.4|) and ( E.5| ). If we assume, as above a translational invariant initial state then the 
truncated equations for the two-point correlation function become, where C(r,t) = < hi(t)h i+r (t) > 

—C(r, t) = -4(1 + 5) C(r, t) + 2( C(r + 1, t) + C(r - 1, t) ) (E.6) 

for r 7^ 1 and 

jC(l, t) = -2(1 + 5 + 7 + a) C(l, t) + 2 C(2, t) (E.7) 

for r = 1. Now equation ( |E.7[ ) can be given the form ( |E.6p if we define C(0,t) = AC(1,£). So now we 



have to solve equation ( E.6|) with this boundary condition. The solution is 

00 

C(r,t) = e- 4 ^ 1 ^[aylr^t) + b y I r+y ^(4t)} (E.8) 

y=l 

where J n (4t) is the modified Bessel function of order n, the a y 's are constants defined by the initial 
distributionF] and b y = Aa y — (1 — A 2 ) J2t=i A y_1_fc a fc . We substitute this in (|E.8| ) and interchange 



16 T 



If S = 2a + 7, eqs. ( E.6| ,E.7) are exact at all times. Then the a y are given by the initial conditions of the full physical 



problem. Otherwise, this relation is considerably more complicated. 
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the sums between y and k. We get after a few manipulations 

oo 

C(r,t) = e- 4 ^ £[a„/ P _,,(4t) + Aa y I r+y ^(4t)} 

y=l 

oo +00 

-(l-A 2 )e- 4 W^ t A- k £ 7^(4*) A*" 1 

k=l j/=— 00 

00 

+ (1-A 2 )e^ 1+S ^Y: a k ]T Iv+r+k-i^t) A^ 1 (E.9) 

fe=l y=—oo 

where we have completed the sum in y in the second term and accordingly subtracted the third term. 
Since e _4t 7 n (4t) < 1 then the sum over y in the third term is convergent if | A | > 1. Hence, provided 
that the sum over a k is convergent this term is bounded by (1 — A 2 ) e~ 4<5t n^Xj k when | A | > 1. 
One can then explicitly perform the sum over y in the second term, using the generating function for 
modified Bessel functions. Hence if | A | > 1 one gets 

00 

C(r,t) = e-^ 1 ^Wr-y^t) + Aaylr+y-tiAt)} 

y=l 



- (l_ A 2) A -r e -l*S-2(A-2 + ^)]t £ ^ 



-k 

a k Ll 
k=l 

00 

+ (1 - A 2 ) e~ A ^ £ a k ^ ■ (E.10) 

k=l y=~oo 



One sees that at long times the first and third terms of this equation decay like e~ 4<5t . Hence if A > 1 
the decay of the two-point correlation functions is dominated by the second term and the long time 
behaviour of the two point correlation function would be proportional to e -[ 45 - 2 ( A ~ 2 + A )]*. Also we see 
that the second term does not appear if | A | < 1, which is in agreement with the result obtained by the 
Bethe Ansatz that there is no bound state with total zero momenta in this range of the parameter A, 
otherwise such state would appear in the decay of correlation functions from a translational invariant 
state. Hence we again obtain the results of Table 2. 

Appendix F. Role of the boundary conditions on the Hamilto- 
nian spectrum 

We discuss the relationship between the quantum Hamiltonian spectra for free and periodic boundary 
conditions. The main point already becomes apparent when just considering asymmetric diffusion 
without any further reactions [^T, [IT], ^7J. While this problem is still very technical for the case of 



a finite particle density, the basic mechanism becomes apparent in a simple way already for a finite 
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number of particles. To make this paper more self-contained, we recall the argument here, following 



For free boundary conditions on a chain of L sites, the quantum Hamiltonian reads 

L-l 



h' = -\ E 



q-q 



where 



A 



q + q 



-i 



1 



q 



(F.l) 



(F.2) 



2 ' 3 V 1 + e 

This is the well-known XXZ quantum chain which is symmetric under the quantum group U q SU(2) 
2^, |36[]. It follows from the Bethe ansatz that for e ^ the energies for L large become El ~ 1 — A -1 



in each sector^]] with iV particles (see also below). We stress that the quantum Hamiltonian of the 
diffusion process on a periodic lattice cannot be obtained by simply taking periodic boundary conditions 
in eq. (FTT). Rather, from the master equation if can be written down directly |2T] 

L 

i 



1 L 



(F.3) 



where a x,y ' z are Pauli matrices and e is related to the diffusion bias via equation ([F.2|) . It is well known 
that the low-lying spectrum of H is gapless (see also below). 

To understand the relation between H and H', we introduce a ± = \(o~ x ± icr y ), and consider the 



non-singular operator 



[/ = exp ( irgYsjVj i ■ 

3=1 



(FA) 



Choosing q = e 2ng , we obtain 

H" = UHU~ l 



2(q + q 



" 3=1 



(F.5) 



which is indeed (almost) the Hamiltonian if'. The distinction comes from the boundary conditions. 
The surface field (q — q^ 1 )(a{ — er£)/8A is absent in if" and one has 



(F.6) 



which make H" non-hermitian, as H already is. It is these unusal boundary conditions which give rise 
to the different properties of H and H', as we now show. 



17 Since H' commutes with the quantum group generators S± and S z , for each state \M) in the M-particle sector with 
energy E' there a (N + M)-particle state \N + M) ~ S N \M) with the same energy, see pi. 
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The calculation of the spectrum of H" proceeds via the Bethe ansatz [59]. First, we consider the 
sector with N = 1 particle. Then the energies are in this sector 

E = 1 - A _1 cos0 , 9 = 2ir(ig + jj (F.7) 

where n is an integer from the set {0, ±1, . . . , +(L/2 — 1), L/2}. For L large, the energies become 

^ cos(2tt(^+ f)) _ /n\ 2 . . ,n 



E = 1 ^— ± ~ 2?r 2 - + 2vri tanh(2vrc/)- + . . . (F.8) 

cosh(27r^) W V ; L v ; 

We see that the real part dtE ~ L -2 and the imaginary part ~ 

Next, we take the sector N = 2. From the Bethe ansatz |6£| we have 



£7 = (2A - cos# - cos 9') A' 1 (F.9) 
0L - 2mgL = 2nl - &(9, 9') , 6'L - 2mgL = 2txV - 0(0', 9) 

where /, I' are distinct half-integers from the set {±|, ±|, . . . , ±^+f-} and 



Asm(\(9-9 ! ) i 

0(9, 9') = 2 arctan T ^ / v (F.10) 



cos (1(6 + 6')) - Acos(§( 
Define 6 = 6 — 2itig, 6' = 6' — 2mg. Then, for small values of the arguments 



(a n/\ 
icoth(2vr^) • ~—~\ ( F - n ) 



which is of order unity. Therefore 



6 = , 9' = ^a' (F.12) 



where a, a' are of order 0(1). It follows that the same cancellation as observed in the sector N = 1 
also takes place here and also that the observed scaling of the energies does not change. Finally, for N 
arbitrary 

E = A- 1 (^NA-J2cos9 n ^ (F.13) 

with 

TV 

9 m L - 2mgL = 2nl m - ]T Q(9 m , 9 n ) (F.14) 

n=l 

and the same argument can be repeated. It follows that for any value of e, the spectrum shows in the 
L — > oo limit no energy gap. 

While this argument shows that for a finite number of particles, the real part of the energies scales as 



L , a much more detailed analysis of the Bethe ansatz equations shows [^T], |27[ that for a finite density, 
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one rather has a scaling $IE ~ L 3 / 2 . This arguments also goes through when particle reactions are 
added P7fl . In the case of more than one species of particles, the relation between H and H' becomes 



even more involved, as explained in 26 



A simple interpretation of this result is as follows. For periodic boundary conditions, the particles 
are free to chase each other. In fact, through a Galilei transformation on can go into the frame where 
the diffusion becomes left-right symmetric, see fTOfl . It is not surprising that it should take a long time 
for the system to reach a steady state (in fact, QE is proportional to the steady state current density) 
and that the relaxation time should grow beyond all bounds with increasing system size L. On the 
other hand, for free boundary conditions the particles cannot leave the system. Rather, they pile up 
on one of the boundaries. The time to reach a steady state is of the order of the time a single particle 
needs to cross the system from one boundary to the other and should be finite. 
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Figure captions 

Fig. 1 Test of the positivity of the reaction rates in the case when the transformation parameters a 
and Y are both real and h 2 > 1. In each subfigure, the shaded regions correspond to the domains 
where one of the following reaction rates is negative with (a) a (b) 7 (c) v and (d) a. 

Fig. 2 Test of the positivity of the reaction rates in the case when a and Y are both real and h = |-\/3. 
Shaded regions correspond to one of the reaction rates being negative with (a) a (b) 7 (c) v and 
(d) a. 

Fig. 3 Test of the positivity of the reaction rates when the transformation parameters A and q are real 
and h 2 > 1. Shaded regions correspond to one of the reaction rates being negative with (a) a (b) 
7 (c) v and (d) a. 
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